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ABSTRACT 

Aims. We investigate the global galaxy evolution over ~12 Gyr (0.05 < z < 4.5), from the far ultraviolet (FUV) luminosity function (LF), 
£^ ' luminosity density (LD), and star formation rate density (SFRD), using the VIMOS-VLT Deep Survey (VVDS), a single deep galaxy redshift 
| survey with a well controlled selection function. 

. Methods. We combine the VVDS Deep (17.5 < I AB < 24.0) and Ultra-Deep (23.00 < i' AB < 24.75) redshift surveys, totalizing ~ 1 1000 galaxies, to 
■ estimate the rest-frame FUV LF and LD, using a wide wavelength range of deep photometry (337 < A < 23 10 nm). We extract the dust attenuation 
of the FUV radiation, embedded in the well-constrained spectral energy distributions. We then derive the dust-corrected SFRD. 
i i ' Results. We find a constant and flat faint-end slope a in the FUV LF at z < 1.7. At z > 1.7, we set a steepening with (1 + z). The absolute 
magnitude M* FUV steadily brightens in the entire range < z < 4.5, and at z > 2 it is on average brighter than in the literature, while <f>* is on 
average smaller. The evolution of our total LD shows a peak atz-2, clearly present also when considering all sources of uncertainty. The SFRD 
history peaks as well atz-2. It first steadily rises by a factor of ~ 6 during 2 Gyr (from z = 4.5 to z = 2), and then decreases by a factor of ~ 12 
| during 10 Gyr down to z = 0.05. This peak is mainly produced by a similar peak within the population of galaxies with -21.5 < M FUV < -19.5. 
Q As times goes by, the total SFRD is dominated by fainter and fainter galaxies. The mean dust attenuation of the global galaxy population rises fast 
I ■ by 1 mag during 2 Gyr from z — 4.5 to z ~ 2, reaches slowly its maximum at z — 1 (A F uv = 2.2 mag), and then decreases by 1 . 1 mag during 7 Gyr 
O ' down to z — 0. 

| Conclusions. We have derived the cosmic SFRD history and the total dust amount in galaxies over a continuous period of ~ 12 Gyr, using a single 
C/3 , homogeneous spectroscopic redshift sample. The presence of a clear peak at z — 2 and a fast rise at z > 2 of the SFRD is compelling for models 
of galaxy formation. This peak is mainly produced by bright galaxies (L > L*_ 2 ), requiring that significant gas reservoirs still exist at this epoch 
' and are probably replenished by cold accretion and wet mergers, while feedback or quenching processes are not yet strong enough to lower the 
_ 1 ' SF. The dust attenuation maximum is reached ~2 Gyr after the SFRD peak, implying a contribution from the intermediate-mass stars to the dust 
, production at z < 2. 

, Key words. Cosmology: observations - Galaxies: evolution - Galaxies: luminosity function, mass function - Galaxies: high-redshift - Galaxies: 
. star formation - ISM: dust 
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0^ ■ 1. Introduction of all processes which shape galaxies along t ime. The study o f 

O ' A , . . , ,. ,„™ TT , the SFH has been pioneered by the CFRS dLillv et all fl996h . 

; A robust determination of the star formation rate history (SFRH) on the evo i u tion of the luminosity density since z ~ 1 , and by 
, is a crucial element to understand galaxy evolution. As time goes Ma< ^^\ (fT^ on the SFR using the CFRS and z ~ 3 - 4 
- ■ by, the reservoirs of pristine gas are transformed into first gener- ^ identified with HSX since ^ ext ensive new SFRD 

<L , ation of stars. During galaxy evolution, feedback processes, gas me asurements have been compiled up to z - 6 (see, for exam- 
^ , accretion, mergers, and environment are expected to affect star p i e ;HopMnI&leacom 200%. 

' formation, and to contribute as well to galaxy mass assembly. ' „ . _ . _ ■ n j • j 

5 . „. nr ,, T ' „ . ' P. J The Star Formation Rate Density (SFRD) is usually derived 

?3 The SFH therefore contains an imprint of the collective outcome J r oo 

from a mean luminosity density, defined as LD = J <p{h) L dL, 



Send offprint requests to: Olga Cucciati (cucciau@oats.inaf.it) with </>(L) being the Luminosity Function (LF) and L a lumi- 

* Based on observations collected at the European Organisation nosity related to the SFR. At first sight, the SFRD is a simple 

for Astronomical Research in the Southern Hemisphere, Chile un- and powerful tool to investigate the cosmic star formation his- 

der programs 072.A-0586 (GTO), 073.A-0647 (GTO) and 177.A-0837 tory. The various data from different samples have persistently 

(LP) at the Very Large Telescope, Paranal, and based on observa- shown a broad picture consistent within factors of ~3, out to 

tions obtained with MegaPrime/MegaCam, a joint project of CFHT high-redshifts (z 6), showing a rise out to z ^ 1 and a decline 

and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) from z ^ ^ whh aR undeaj . evolution within 1 < z < 2 .5 be- 

which is operated by the National Research Council (NRC) of Canada, , , . , , ..^ . , — 

. . , , „ . , „ TT . ' XT ., I cause this redshift desert has remained difficult to probe. From 

the Institut National des Science de 1 Umvers of the Centre National , , , , . _ , o^t,,^ , , 

de la Recherche Scientifique (CNRS) of France, and the University the P resent day to z - 1, a steady rise of the SFRD by one order 

of Hawaii. This work is based in part on data products produced of magnitude is firmly corroborated using various calibrators of 

at TERAPIX and the Canadian Astronomy Data Centre as part of SFR - like far ultraviolet (FUV), far infrared (FIR), Ha, radio 

the Canada-France-Hawaii Telescope Legacy Survey, a collaborative - but the scatter amongst different measurements remains large. 

project of NRC and CNRS. Because of the number of uncertainties that remain along the 
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chain of transformations to be applied to galaxy counts and lu- 
minosities to be converted into star formation rates, the exact 
shape of the SFRD evolution still remains to be established. The 
selection function of each galaxy sample, including the imaging 
surveys depth and image quality or the redshift completeness, re- 
quires a number of corrections to compute volume densities, and 
the complete shape of the luminosity function remains specula- 
tive, particularly at the faint end. Moreover, the transformation of 
luminosities to SFR depends on assumptions on the conversion 
factors and on the amount of dust attenuation, still open issues. 

To be able to derive a coherent galaxy evolution model of 
galaxy evolution, it is necessary to trace the SFRD evolution 
with the same reference indicator throughout cosmic time, and 
within a deep galaxy redshift survey with a simple, well defined 
and controlled selection function. 

The selection of spectroscopic targets may have a strong im- 
pact on the type of galaxies studied and hence on the relevance of 
this sample to the full galaxy population. A tight control of the 
selection function is necessary to avoid the propagation of bi- 
ases to some types of galaxies or redshift ranges. Several meth- 
ods are used to pre-select galaxies, either using colour criteria 
tailored to different redshift ranges and populations (e.g. LBG, 
BzK, radio, far-IR,...), or pure magnitude selection. A follow-up 
spectroscopic survey with high success rate has then the advan- 
tage to provide a sample of galaxies with known redshift and 
controlled uncertainties, with stars and broad-line active galaxy 
nuclei clearly identified as pre-selection techniques are largely 
unable to fully discard them. 

Using rest-frame ultraviolet galaxy luminosities has become 
a common approach at high redshift, as it is applicable up to th e 
highest redshifts studied so far (z - 7, e.g. lBouwens et al l2009h . 
The non-ionising ultraviolet light (912-3000 A) is emitted by 
relatively massive (> 3 M Q ), short-lived (< 3 10 8 yr) stars and 
it traces the SFR averaged over the last ~ 10 8 yr once it is cor- 
rected for dust attenuation dKennicuttll998h . This SFR estimator 
assumes a constant SFR over longer times than the very massive 
stellar population (>> 15 M Q ) contributing to the Ha emission. 
As opposed to the Ho- instantaneous SFR estimator, it does not 
disentangle whether the radiation is linked to the creation of new 
stars (< 10 7 yr) or to more evolved stars (< 10 8 yr), and thus i t 
is sensitive to ageing of star formation regions (ICalzettill200i . 
An ideal measurement of the instantaneous SFR, not affected by 
uncertainties on dust, would be based on the simultaneous use 
of Ha and infrared emissivities, since the ab sorbed ionising flux 
heats the dust which re-emits in the infrared ( ICalzettil2 008). The 
difficulty is that it has not yet been possible to assemble large, 
deep and statistically complete samples for both estimators. 

Here we aim to derive the SFRD evolution since z ~ 4.5 us- 
ing the rest-frame UV luminosity density from the VIMOS VLT 
Deep Survey (VVDS) samples 'Deep' and 'Ultra-Deep', total- 
izing ~ 11000 galaxies with spectroscopic redshifts. We take 
advantage of being able to compute rest-frame ultraviolet emis- 
sivities over a large and deep area of sky, which enables to trace 
its evolution in a consistent way since z — 4.5. 

This paper is organised as follows. Section 2 gives a sum- 
mary of our VVDS data sets. Sections 3 & 4 describe the 
ultraviolet luminosity functions (LF) and densities (LD) from 
z = 0.05 to z — 4.5 as derived from these data. Section 5 
presents our results for the evolution of the dust attenuation and 
the dust-corrected SFRD in the past ~ 12 Gyr. Sect. 6 sum- 
marises and discusses our results. Technical details are given 
in the Appendixes. Throughout this paper, we use the AB flux 
normalisation, and we adopt the concordance cosmology (Cl m , 



h) = (0.3, 0.7, 0.7), with which the age of the Universe is 
5.7,3.2,1.3 Gyr at z= 1,2,4.5. 

2. Data 

2.1. The VVDS Deep & Ultra-Deep spectroscopic surveys 

The VIMOS-VLT Deep Survey (VVDS) is a spectroscopic 
investigation of distant sources, carried out with the high- 
multiplex, wide-field (224 arcmin 2 ) Visib le Multi-Object 
Spectrograph (VIMOS. ILe Fevre etakl I2003I) mounted at the 
Nasmyth focus of MELIPAL, the third of the four 8.2m ESO- 
VLT Unit Telescopes in Paranal, Chile. The VVDS is composed 
of three /-selected surveys totalizing about 47000 spectra of 
gala xies, quasars, and stars; (1) a Wide survey (17.5 < Iab ^ 
22.5 iGarilli et al.ll2008h, (2) a Deep survey (17.5 < I AB < 24.0 
ILe Fevre et al.l l2004bl l2005bl) . and (3) an Ultra-Deep survey 
(23.00 < i' AB < 24.75, Le Fevre et al., in prep.). It spans a wide 
redshift range; from z > to z — 5 for tar geted sources, and up 
to z — 6.6 for serendipitous Lya emitters ( ICassata et al.ll2.01 lb . 
The VVDS surveys are purely flux limited, and are free from 
any colour pre-selection or galaxy-star separation when prepar- 
ing the target catalogues. In the Deep survey, the projection of 
galaxy size on the x-axis of the image has been used to maximise 



the number of targets in each VIMOS pointin g (see lllbert et al 
2005 for further details). We refer the reader to lMcCracken et al 



(120031) for a complete discussion on photometric completeness 
as a function of Iab magnitude and surface brightness. The con- 
clusion of their analysis is that the VIMOS Deep Imaging Survey 
"is essentially free of selection effects until at least Iab = 25". 
Namely, considering the flux limit of the Ultra-Deep survey, 
Iab ^ 24.75, our photometric catalogue is 90% and 70% com- 
plete in surface brightness down to 24.5 and 25 mag/arcsec 2 , 
respectively. 

In this study, we use both the Deep and Ultra-Deep surveys 
obtained in the VVDS-0226-04 field. The Deep spectra were col- 
lected with integration times of 16000 seconds using the LR-Red 
grism (R = 210, 5500 < A < 9500 A), over 2200 arcmin 2 of sky 
area. Th ey include the 9842 spectra described in lLe Fevre et al.l 
(l2005bl) . plus 2826 spectra acquired later with the same set-up 
(Le Fevre et al. in prep.). With the Ultra-Deep survey we ob- 
tained repeated observations of several z > 1 .4 targets from the 
Deep survey with assigned VVDS quality flags of 0, 1, or 2 (i.e., 
with the lowest confidence level in the spectroscopic identifi- 
cation) in order to assess their real redshift distribution as de- 
tailed in Le Fevre et al. (in prep.). The Ultra-Deep spectra were 
collected using both the LR-Red grism and the LR-Blue grism 
(3700 < A < 6700 A, R = 180) with integration times of 
65000 seconds in each set-up. They consists of 1200 new targets 
acquired over 576 arcmin 2 of sky area within the 2200 arcmin 2 
sampled by the Deep survey. Our present work is based on 10141 
and 622 galaxy spectra at 0.05 < z < 4.50 from the Deep and 
Ultra-Deep surveys, respectively. Those spectra have a spectro- 
scopic identification at a confidence level higher than ~ 50%, 
60%, 81%, 97% and 99% (corresponding to the VVDS quality 
flags 1,9, 2, 3 and 4). 

We do not have a measured redshift for every source to a 
given apparent magnitude limit in the observed field of view; 
the averaged target sampling rates are ~ 24% and ~ 4% for the 
Deep (17.5 < Iab — 24.0) and Ultra-Deep (23.00 < i' AB < 24.75) 
surveys, respectively. For each survey, we accurately estimated 
the selection function accounting for the facts that (i) a fraction 
of sources of the parent photometric catalogue was targeted for 
spectroscopic observations; (ii) a fraction of the targeted sources 
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yields reliable redshifts. The computation of the selection func- 
tions practically translates into weights applied to each galaxy 
when computing the luminosity functions. In addition, in the 
Deep survey, we used the results from our repeated observations, 
to obtain a 100% confidence level in the redshift measurement 
of lower quality spectra. The weighting schemes for the Deep 
and Ultra-Deep surveys are independent, and we refer the reader 
to Appendix lA. 1 l and Appendix IA.2l for a detailed description. 

2.2. The absolute magnitudes 

For the astrometric and photometric catalogues, we take the 
following broadband imaging surveys over the VVDS-0226-04 
field acquired at the Canada-France-Hawaii Telescope (CFHT) 
using w ide-field mosaic camer a s: the BVI VVDS im aging 
survey dMcCracken et alj 120031: |Ee Fevre et al. | l2004ah with 
the CFHT-12K camera, the u*g'r'i'z' CFHT-Legacy Survey 
(CFHTLS-D1/W1 field, T0005 releas e) with the MEGA CAM 
camera and the JHK S WIRDS survey dBielbv et alj|201 lb with 
the WIRCAM camera. 

To derive intrinsic lumin osities, we use ou r Algorithm for 
Luminosity Function (ALF, Ilbert et aT| 120051) . that integrates 
routines of the code Le Phar$\. Absolute magnitude measure- 
ments are optimised accounting for the full information given 
by the above multi-band photometric data, in a way which min- 
imises the dependency o n the templates c hosen to fit the ob- 
served colours (se e A. 1 in Ilbert et al. |2005|). W e choose the tem- 
plate library from iBruzual & CharlotM 2003b modulated by the 
attenuation of the intrinsic stellar continuum, A(A) = k(A)E(B - 
V). We consider a grid for the intrinsic colour excess of the 
continuum light (E{B - V)=0, 0.1, 0.2, 0.3, 0.4 and 0.5 mag), 
an d an empirical attenu ation curve, k(A), of the form defined 
by ICalzetti et al.l d2000t) . We verified that the absolute magni- 
tudes are indeed mostly independent from the chosen template 
set, especially for those derived in the rest-frame FUV, NUV, 
U and B photometric bands. For this purpose, we computed 
them a gain using different template sets , like PEGASE.2 tem- 
plates {Roc & Rocca-Volmerangelll997b or a mix ed set with 
IBruzual & Chariot! d2003l) a nd iPolletta et al.l (120071) templates, 
as used in lllbert et al.l (l2009h . This very weak dependency of the 
derived absolute magnitudes on the adopted templates is mainly 
due to the fact that we use i) a wide (337 < A < 2310 nm) 
wavelength range with the observed u*Bg'Vr'i'Iz'JHK s photo- 
metric broadbands to obtain a robust template fit, and if) the u*- 
CFHTLS broadband filter (our bluest photometric information, 
337 < A < 41 1 nm) that begins to sample the non-ionising ultra- 
violet continuum (9 1 .2 - 300 nm) at z > 0. 1 . The u* photometric 
information is deep enough to be complete down to our VVDS 
spectroscopic limiting magnitudes, based on the primary selec- 
tion from the 7-band"cFHT-12K images (725 < A < 930 nm). 
This is the case at least up to z — 3.5 when the Lyman-break 
feature is shifted towards our redder photometric bands, i.e., 
Bg'Vr'i'I. The observed m*-381 corresponds to the rest-frame 
NUV-250 at z 0.4 and to the rest-frame FUV-150 at z 1.5. It 
spans the NUV-250 band (184 < A < 280 nm) from z = 0.2 to 
Z= 1.2 and the FUV-150 band (135 <A < 175 nm) from z = 0.9 
to z = 2. At z > 2, the rest-frame FUV-150 and NUV-250 lu- 
minosities are covered by observations with filters redder than 
m* (up to the J band at z = 4 for the NUV-250). Our rest-frame 
FUV absolute magnitudes are based on template extrapolation 
only for z < 1, but still we verified that our FUV-based results 
agree with our NUV-based ones. As an example of these consis- 



tency tests, in Appendix|B]we compare the rest-frame NUV and 
FUV luminosity functions. 



3. The ultraviolet Deep & Ultra-Deep luminosity 
functions 

In this work, we have primarily estimated the rest-frame FUV 
luminosity functions and densities to derive the SFRD history. 
The recent SFR is traced by the intrinsic non-ionising ultravio- 
let stellar continuum (91.2 - 300 nm) of galaxies (see Sect.[TJ. 
Within this UV range, the far UV radiation (FUV-150) is a better 
SFR indicator than the near UV radiation (NUV-250), because 
the NUV is contaminated by evolved stars, while the FUV is 
dominate d by the radiation f rom new, massive, short-lived stars 
(see, e.g jMadauetal .1119981) . 

3. 1 . The method 
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The ga laxy luminosity function (LF) usually follows a lSchechterl 
d 1976b function characterised by a break luminosity, L*, a faint- 
end slope, a, and a normalisation density parameter, </>*. The 
LF is a fundamental measurement of the statistical properties 
of the galaxy population; it is the distribution of the galaxy co- 
moving number density as a function of their intrinsic lumi- 
nosity at a given epoch. Despite its simple definition, its esti- 
mation requires careful analyses of the survey strategy, the se- 
lection criteria, and the completeness. We derive it using our 
code ALF that includes the non-parametric 1 /V max , SWML, and 
C + and the par ametric STY lumi nosity function estimators (see 
Appendixes in lllbert et alj|2005l and references therein). Each 
estimator presents advantages and drawbacks, and each one is 
affected differently by different visibility limits for the various 
galaxy types detected in deep flux-limited surveys. In a given 
observed band, galaxies are not equally visible to the same abso- 
lute magnitude limit mainly due to the spectral type dependency 
of the K corrections (i.e., different spectral energy distribution 
-SED- for different spectral type galaxies). Within a given red- 
shift range, the use of several non-parametric LF estimators al- 
lows us to empirically determine the absolute magnitude range 
in which galaxies are equally visible within our spectroscopic 
surveys. The luminosity limit down to which all g alaxy popu- 
lations are visible is called LF bias. As discussed in lllbert et al.l 
(l2004h . 1/V max and C + methods are affected by this bias in a 
different way than the SWML and STY methods. If 1 /V max or 
C + at a given luminosity starts giving different results from the 
SWML or STY, it means that a bias in the global LF is present. 
The bias is significant if the two estimators differ of more than 
the statistical uncertainties (Poisson errors, in our case). In each 
redshift bin explored, we set the LF bias at the brightest abso- 
lute magnitude where l/V maA - or C + are different by more than 
1-cr from SWML or STY. Our LF parameters are estimated with 
data brighter than the LF bias limit, to derive an unbiased LF 
faint-end slope of the global population. This is particularly im- 
portant when com bining surveys of different depth. We use the 
ISchechterld 1976b functional form for the STY estimate to calcu- 
late the Schechter parameters, because the results are more ro- 
bust than those with a simple fit with a Schechter function of the 
non-parametric LF data points. The resulting faint-end slope, a, 
is independent on the luminosity binning, and since the LF pa- 
rameters are highly correlated to each other, we can account for 
the allowed range for each Schechter parameter as derived by 
the likelihood, in addition to the Poisson uncertainties typically 
quoted. 
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Fig. 1. Rest-frame FUV-band luminosity functions in ten redshift bins from z = 0.05 to z — 4.50 of the merged Deep+Ultra-Deep 
VVDS sample. Redshift ranges are indicated in each panel. Red circles represent the 1 / V ma x data points with Poisson error bars, and 
the number of galaxies is given in parenthesis. 1 /V max points are plotted up to the LF bias limit (see Sect. [XT) of the merged sample, 
represented with a vertical blue short-dashed line. This line corresponds also to the LF bias limit of the Ultra-Deep sample. The 
vertical cyan long-dashed line corresponds to the LF bias limit of the Deep survey. There is no Ultra-Deep data at 0.05 < z < 0.2, 
so the only bias limit is the one of the Deep sample. The black solid curve and its associated orange shaded area is the STY LF 
estimate assuming free Schechter parameters and its associated error. At 3.5 < z < 4.5, the STY LF fit does not converge using all 
free parameters. In the first redshift bin, the dashed curve is the STY LF estimate when setting M* FUV = -18.12 (see text), and for 
reference it is reported as a dotted curve in the other panels. In the 1.7 < z < 2.5, 2.5 < z < 3.5, 3.5 < z < 4.5 panels, the dashed 
curve is the STY LF estimate when setting a faint-end slope that evolves with increasing redshift (i.e., a = -1.3, -1.5, -1.73, 
respectively). The Schechter parameters of all the curves shown in this plot are listed in Table IB. fl while Table [1] summarises those 
that we retained for our study. 



3.2. The UV luminosity functions 

We compute the rest-frame FUV LFs using a unique merged cat- 
alogue which includes both the Deep and Ultra-Deep surveys, 
for a total covered magnitude range of 17.5 < Iab ^ 24.75. This 
way, we exploit both the large magnitude range covered by the 
Deep survey and the depth reached by the Ultra-Deep survey. 
This leads us to a robust determination of the LF shape and nor- 
malisation. 

When using flux-limited surveys with various apparent lumi- 
nosity depths, a coherent weighting scheme must take into ac- 
count the possible overlap of the flux ranges covered by the dif- 



ferent surveys. We detail the weights to be applied in this case in 
Appendix lA.3l and we verify their robustness in Appendix I A.4I 
Fig. Q] shows the rest-frame FUV-band LFs (l/V max and STY) 
obtained with the merged catalogue with 17.50 < Iab ^ 24.75, 
in ten independent redshift bins from z = 0.05 to z = 4.5. The 
non-parametric 1 1 Vmax data points are plotted up to the LF bias 
limit, that is where those from the SWML and C + methods are 
in agreement (the latter estimates are not shown throughout the 
paper for clarity in the figures). At z < 0.2, due to the size of our 
survey fields, rare bright nearby galaxies are not observed and 
the bright-end of the LF cannot be constrained as shown by the 
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Table 1. LF Schechter parameters, LF absolute magnitude bias limits, LDs, SFRDs and dust attenuation derived from our rest-frame 
FUV-150 LF estimates in 10 redshift bins from z = 0.05 to z = 4.5 assuming (Q„„ Q A , h) = (0.3, 0.7, 0.7). 
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Notes. Values derived for the merged VVDS Deep+Ultra-Deep sample at 17.50 < I AB < 24.75. There is no Ultra-Deep data at z < 0.2. (D) Redshift 
bins. Schechter parameters (M* , a and (f>*) of our VVDS rest-frame FUV-150 luminosity function that we use throughout this work. (c) Absolute 
magnitude limit in FUV down to which the survey is complete is terms of galaxy types (see Sect. l3.lt . for the Deep and the Ultra-Deep surveys, 
respectively. W) Rest-frame FUV luminosity density uncorrected for dust. The quoted uncertainty includes errors on the LD induced by the STY 
LF fit, the Poisson noise, the cosmic variance, and the weighting scheme. When a is set at z > 1.7, the STY errors include also the percentage 
uncertainty accounting for the two extreme values of a (-1.1, -1.73) with respect to those quoted in the third column (see text). The relative 
contribution of each source of error is detailed in Table [2]. {e) Star formation rate density uncorrected for dust, derived from the rest-frame FUV 
luminosity density. For the errors see note (d). lt) Average dust attenuation in FUV. ' s > Star formation rate density corrected for the dust attenuation 
A FUV . For the errors see note (d). 



large error (shaded area in Fig.Q]) associated to the STY estimate. 
For this reason, in the first r edshift bin we set M* to the local 
FUV value, M* FUV = -18.12 dWvder et al.ll2005l) . Conversely, at 
z > 1.7, the faint-end of the LF starts to be loosely constrained. 
Recent results in the literature (see Table IB.2b seem to indicate 
a steep faint-end slope of the rest-frame FUV-band LF at z > 2. 
We therefore opt for a slope set to evolve with increasing red- 
shift, using the p arameterisation, a(z) = A(l + z) + B (see, e.g., 
iRvan et al]|2007l) . To derive A and B, we use our highest red- 
shift estimate of a (a = -1.1 at z ~ 1.45, see Table IB. It and 
a = -1.73 at z = 3.8 taken from lBouwens et alJ d2007l) . In this 
way, we obtain a = -1.3, -1.5, and -1.73 at 1.7 < z < 2.5, 
2.5 < z < 3.5, and 3.5 < z < 4.5, respectively. In addition, in 
all the redshift bins at 1.7 < z < 4.5 we compute <p* and M* 
using also extreme non-evolving a cases: a = -1.1 as estimated 
at 1.2 < z < 17, an d a - -1.73, as estimated at z ~ 3.8 by 
iBouwens et al.l d2007l) . Table IB. ll summarises the Schechter pa- 
rameters for all above cases. 

From now on, we consider our final best Schechter param- 
eters ((p*, M* , a) of the rest-frame FUV LFs as those esti- 
mated with the STY method with M* set at the local value at 
0.05 < z < 0.2, with all free parameters at 0.2 < z < 1.7, 
and with a set to evolve at 1.7 < z < 4.5. They are reported 
in Table [1] and the corresponding LFs are plotted all together in 
Fig. [2] They are used to derive the luminosity densities (Sect. |4]i 
and the star formation rate densities (Sect. |5j. At z > 1.7, the 
uncertainties of these densities will include the possible span of 
a between the values -1.1 and -1.73. 



3.3. Main features of the derived FUV LF 

Our results show a constant and flat (a -1) FUV faint-end 
slope at z < 2, while we set it to steepen with z at z > 2, where 
our data can not constrain it. We find that M* brightens mono- 
tonically with increasing redshift for the entire redshift range 
explored. Conversely, at z > 0.9 (f>* starts decreasing, with a re- 
markable drop especially at z > 2. These trends are evident in 
Fig.|2]and in Table[T] and are also summarised in Figj8] 



Table 2. Relative uncertainty to the total LD (sixth column of 
TableQ]) from different sources of error, expressed in percentage. 
Their sum in quadrature give the total error quoted in TableQ] 



Az" 



0.05-0. 
0.2-0 
0.4-0 
0.6-0 
0.8-1 
1.0-1 
1.2-1 
1.7-2 
2.5-3 
3.5-4 



STY' 



a range c Po 
errors in % 



CV 



+ 10.9 

-<->.» 
+7.9 



+53 
-4.9 
+3.5 
-3.3 
+7.3 

-6 4 





6.3 


40 


~0 




3.0 


26 


~0 




2.9 


22 


~0 




2.6 


20 


~0 




2.8 


19 


~0 




3.4 


19 


~0 




3.4 


12 


~0 


+ 112 
-12 

+56 

-32 

+0 

-55 


7.2 
5.8 
8.8 


10 
9 
10 


+ 11.3 
-7.8 
+ 17.0 

-46 9 



Notes. For details on the computation of each error see Sec. 14.11 
<n) Redshift bins. " J) Relative error induced by the STY LF fit. (c) Relative 
uncertainty derived from the difference in the total LD between the 
'best' a value quoted in Table [T] and the lower (a = -1.1) and up- 
per (a = -1.73) limits considered. It is considered only when a is set. 
(<,) Relative error induced by Poisson noise. w Relative error induced by 
cosmic variance. ( ^ Relative error derived accounting for the error of 
the galaxy weighting scheme in the LF fit. 



3.3.1 . A persistent flat faint-end slope at z < 2 

Table IB~2l lists a values used in the literature. Contrarily to our 
work, at z < 2 most previous studies were not faint enough 
to determine the faint-end slope, which was then often set to 
values as different as —1.1 < a < -1.6. Neverthele ss, for in- 
stance, lArnouts et all d2005l) and lOesch et ail d2010l) estimated 
a steep, but also constant, FUV faint-end slope a =* -1.6 at 
0.4 < z < 1.2 and a = -1.7 at 0.75 < z < 2.5, respectively. 
In co ntrast, low-z values ( z < 0.4) are found r ather flat a - -1.2 
(e.g.. lArnouts et all 120051 IWvder et al.ll2005l) . To establish the 
robustness of our results, Appendix|B]details our three tests sum- 
marised below. 
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Fig. 2. The VVDS rest-frame FUV-band luminosity functions 
from z - 0.1 up to z - 4 (from the thinnest to the thickest 
curve), fitted with a Schechter functional form parametrised with 
the Schechter parameters (a, </>*, M*) reported in Table Q] Lines 
are dashed for magnitudes fainter than the LF magnitude bias in 
each redshift bin (see Sec. 13. Q . 

For the first test, we derived the rest-frame FUV 1 /V max esti- 
mates using the VVDS photometric redshift catalogue, 1 .25 mag 
deeper than the VVDS spectroscopic catalogue (i.e., Iab = 26), 
over the same VVDS area for a total of ~ 43000 sources (see 
Fig. IB. lb . A similar flat faint-end slope is found. Thus our de- 
rived flat slopes are not due to a large amount of faint blue 
galaxies missing in the spectroscopic sample, or to an inadequate 
weighting scheme. With this multi-wavelength {u*g'r'i'z'JHK. s ) 
catalogue, we also tested the effect of selecting galaxies upon 
their observed M*-band up to 26 mag (i.e, based on a rest-frame 
ultraviolet selection at z < 1) rather than the VVDS observed 
/-band selection (i.e., based on a rest-frame optical selection at 
z < 1). Even though the faint-end slope of such a sample is 
slightly steeper, but not by more than 0.1, it stays very constant 
and close to flatness at z S= 1 . Thus, our /-band selection is un- 
likely missing a significant population of rest-frame blue galax- 
ies, that could have caused an hypothetic steeper slope. 

In the second test, we confirm that the blue, or very blue 
galaxy population of the VVDS /-selected spectroscopic survey, 
does not give a steeper slope than the global galaxy population 
(see Fig. IB. 2b . In rest-frame bands redder than UV (e.g. U, B), 
the shape of the LF i s highly dep e ndent upon the type of galax- 
ies. For instance, in IZucca et aD d2006l) at 0.4 < z < 0.9, we 
found a - — 0.9(/-l) in t/(/Z?)-band for the blue galaxy popu- 
lation, a =s — 1.7(/-1.6) in t/(//?)-band for the very blue galaxy 
population, and a —1.2 in both U and B bands for the global 
population. The diversity of the shape occurs because in bands 
redder than FUV, the radiation is contaminated by a long-lived 
stellar population, that makes highly distinctive the galaxy spec- 
tral energy distributions for different galaxy types. In contrast, 
the FUV-band radiation is dominated by a single stellar popula- 
tion of short-lived, massive stars, whatever the type of galaxies. 
This explains why we observe little change in our a FUV slope 
between the global galaxy population and the blue, or very blue, 
galaxy populations. Another possible additional effect is the ris- 



ing d ust attenuation towards shorter wavelengths dDriver et al.l 
120081) . that increases the dispersion of the rest-frame B - FUV 
colours. This dispersion, if dependent on the luminosity, may 
weaken the B-band steep slope when using the same galaxy sam- 
ple to compute the LF in FC/V-band. 

In the third test, we derive the rest-frame NUV LFs at 0.05 < 
z < 4.5 (see Fig. lB.3l and Table IB. ft . The faint-end slope is also 
flat at z ^ 1 .5 also in the NUV band, showing that the flatness of 
our FUV slope is not due to possible problems induced by our 
template extrapolations at z < 1 . 

In summary, the robustness of our flat UV faint-end slope 
is well-established at z < 2, and it is more consis tent with the 
faint- end slope the local universe, i.e. a -1.2 (IWvder et al.l 
120051) . Despite its low value, we will see (Sec. 15.2b that our dust- 
corrected SFRD using our rest-frame FUV LFs fully agrees with 
other measurements at z < 1, a cosmic period where the SFRD 
seems now well-constrained with different SFR estimators. This 
is because we correct the FUV radiation for dust attenuation in 
an appropriate way as we describe later. 

3.3.2. Comparing FUV LFs at z > 2 

At z > 2, M* FUV keeps on brightening, while <f>* drops. Thanks 
to the depth of our data, these trends are robust, whichever the 
choice of a (see Table IB. ll . 

We remark that the drop in <f>* that we see at z > 2 is likely 
unrelated to incompleteness of low surface brightness objects, 
although this effect can not be completely ruled out. As men- 
tioned in Sec. 12.11 the VVDS survey photometric catalogue, 
down to Iab ~ 24.75, is 90(/70)% complete in observed sur- 
face brightness down to 24.5(/25) mag/arcsec 2 . As discussed in 
iMcCracken et alJ (120031) (Sec. 4. 1), objects with particularly low 
surface brightness, that could fall below our detection limit at 
high redshift, are very few, while 'normal' galaxies would be in 
any case detected. Still, a possible concern can be raised if we 
consider that the above-mentioned observed surface brightness 
limits correspond to brighter and brighter absolute limits going 
to high redshift. We will see (Sec. I5.3l > that the major contribu- 
tion to the total SFRD at z > 3 is given by very luminous galax- 
ies, so in principle our results at these epochs should be driven 
by galaxies with an absolute surface brightness bright enough 
to be within the detection limits. Nevertheless, we could miss a 
population of low absolute surface brightness galaxies, for which 
we would not be complete at high redshift. We performed a very 
simple test to quantify, at a first order of magnitude, this possible 
effect. We considered the redshift bin 3.5 < z < 4.5, where the 
I-band filter, used to select VVDS galaxies, corresponds to the 
rest-frame FUV light. T his way, to pass from the surface bright- 
ness in I-band shown in IMcCracken et all (12003) to an absolute 
FUV-band surface brightness we can neglect K— and colour- 
corrections. We verified that all the galaxies brighter than the 
DeepQ LF bias (Mfuv - -22) have an absolute surface bright- 
ness brighter than 18 mag/arcsec 2 , that coresponds, at z ~ 4, to 
the observed limit of 25 mag/arcsec 2 (considering only the cos- 
mological dimming, with no K— or colour-corrections), and only 
a small fraction is fainter than 17.5 mag/arcsec 2 , corresponding 
at z ~ 4 to the observed limit of 24.5 mag/arcsec 2 . In summary, 
even if small effects can not be completely ruled out, we believe 
that our results are robust with respect to surface brightness com- 
pleteness. 



2 In this redshift range, only a very small fraction of the galaxies that 
we use to compute the LF are from the Ultra-Deep sample, as can be 
seen in Fig. IA.2l 
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We also verified that the low <p* value at z > 2 is robust with 
respect to the weighting scheme (see Appendix[A} that we apply 
to the spectroscopic sample. Thanks to this weighting scheme, 
we are indeed able to properly recover the total galaxy popula- 
tion in the photometric catalogue. Fig. lB.ll shows the FUV-band 
LFs computed using the entire photometric catalogue (with pho- 
tometric redshifts): it overlaps with the LFs computed using the 
VVDS Deep+Ultra-Deep spectroscopic sample, reinforcing the 
reliability of our weighting scheme. 

Table IB. 21 shows clearly that the determination of a is a ma- 
jor source of uncertainty, and even when it is fitted and not set, 
it can be as flat as ~ —1.2 or as steep as ~ —1.7. The scatter 
on a is accompanied by analogous ones for M* FUV and <p*. We 
verified that our M* FUV is generally brighter than other works in 
literature, at all z = 2,3,4, wh ile our 0* is much sm aller, espe- 
cially at z — 3 and 4 (see e.g. IGabasch et al.ll2.004l a t z = 2, 3, 
Steidel et alj|1999l a t z = 3. 4, lArnouts et al.1 12005) at z = 3, 



Bouwens et al.ll2007l at z = 4 iReddv & Steidell(2009( at z = 2. 3. 
Paltani et al.ll2007l at z = 3 .5 IGabasch et al.ll2004lat z = 3.5, 4 5, 
Sawicki & Thompsonll2006l at z. = 2. 3, 4, iTresse et alll2007l at 



z = 3, 3.6,4.3). We remark, however, that in a few cases we are 
in agreement with some previous work for what c oncerns M* FUV 
and/or </>* : we agree on both M* FVV and 0* with IGabasch et al.l 
d2004 at z = 2.5 and w ithlOesch et alJl2010l at z ~ 2, and 
we agree on M' F[ ,„ with iPaltani et al.1 l2007T at z — 3.5 and 
lArnouts et aT1l2005l at z ~ 3. 

This comparison cannot be exhaustive, because it is hard to 
directly compare the LF shapes from different works. This hap- 
pens for the following reasons. First, the three Schechter param- 
eters are linked together; therefore, setting a will in some sense 
determine also M* and <p*, and we already showed that a have 
been set, or determined, to very different values. Second, com- 
puting the Schechter parameters fitting the 1 / V mflA values or us- 
ing the STY method can lead to slightly different results. Third, 
different FUV-band filters or central wavelength have been used 
in the literature (peaked at 1350 < A < 1700, see TablelB.2b. 

As a final remark about comparing LFs in the literature, we 
detail in Appendix ICl the difficulty in comparing FUV luminos- 
ity functions directly-obse rved (as ours) with those derived from 
FUV number counts (as in lSteidel et al.l 19991 for instance). With 
this analysis, we explain t he following remarkab le evidence: at 
2.7 < z < 3.4, on one side iLe Fevre et al" and Le Fevre 

et al. (in prep.) show that the VVDS presents surface number 
counts (per ap parent magnitude) at least two times larger than 
those found by ISteidel et al.l (1 19991) . and this is particularly ev- 
ident for bright galaxies (Lib < 23.5); on the other side, the 
VVDS rest-frame FUV LF is only ~ 50% higher at bright mag- 
nitudes than the one found bv ISteidel et ail d 19991) . We refer the 
reader to Appendix ICl for details. The conclusion of this exer- 
cise is that number counts, in colour-selected samples, are less 
representative of the complete galaxy population than those in 
flux-limited samples, as the respective n{z) of these samples are 
very different. This exercise also confirms that one can not easily 
transform number counts to LF in the case of a skewed n(z). 



4. The FUV comoving luminosity densities 

We derive the mean comoving luminosity density (LD) in each 
redshift bin as LD = J h bngh ' cf>(L) L dL, where <p(L) is the lu- 
minosity function assuming a ISchechterl(fl976l) functional form 
as done in Sect. [3] We set L/ m „, = 10 15 W Hz 1 and Lh r i g ht = 
10 25 W Hz 1 (corresponding to Mf a - m , ~ -3.4 and Mhrieht ~ 
-28.4), to adopt the same limits as in the compilation o flHopkinsI 



(12004 and to allow an easier comparison with the literature. 
Given the typical shape of the LF (steeply decreasing for bright 
galaxies), setting Lbnght as above or Lt, r i g i lt = °° does not make 
any difference. Setting the faint limit at Lfaint — 

10 15 W Hz" 1 

gives a total LD on average 0.3% different from the definition 
With Lf a i m = 0. 

4. 1 . The VVDS FUV comoving luminosity density 

The left panel of Fig.|3]shows the FUV-band LD that we derived 
from our best LF determinations at < z < 4.5 and tabulated 
with their total uncertainties in TableQ] Our LD uncertainties in- 
clude errors from the STY LF fit, errors due to cosmic variance, 
Poisson noise, and errors associated to our weighting scheme. 
The relative uncertainties induced on the total LD by each dif- 
ferent source of error are listed in Table|2] 

Errors derived from the STY estimator are underestimated at 
z > 1.7 since a has been fixed in these redshift bins. In these 
cases, we computed the LD also assuming the two extreme val- 
ues of a (-1.1,-1.73), and the percentage difference with our best 
LDs at z > 1.7 was added in quadrature to the STY errors. 
Cosmic variance (CV) erro rs are computed using the recipe in 
iDriver & Robothaml (20101). We com puted CV errors also with 
the recipe in Trenti & Stiavellil (120081) . and we found lower esti- 
mates at z < 1 .2, going from ~ 20% at z — 0. 1 down to ~ 12% 
at z - 1.1. The same computation based on Somerville et al.1 



(2004) gives results similar to those bv lTrenti & Stiavellil d2008l 
We decid ed to use the more conserva tive values found with the 
recipe in IDriver & Robothaml (12.0 1 Oh . Errors associated to the 
weighting scheme are computed as follows. The weight for each 
galaxy has an error derived in propagating the Poisson noise in 
the weighting formulas detailed in Appendix lAl We computed 
again the LFs using the weights plus and minus this error, and 
we derived the maximal and minimal LDs. Their percentage dif- 
ference with respect to the best LDs gives errors arising from the 
weights. This error is negligible below z ~ 1.7. In each redshift 
bin, the final LD uncertainty is the addition in quadrature of all 
these sources of errors. From Table |2l it is evident that CV er- 
rors dominate at z < 1 .7, while at z > 1.7 the total uncertainty 
is dominated by errors from the weighting scheme and by the 
percentage difference when assuming the two extreme values of 
a with respect to the chosen a value. It is worth noticing that 
our largest uncertainty, the upper error bar at z ~ 2, is given 
mainly by the uncertainty on the value of a: it includes the dif- 
ference between the LD with our best a (= -1.3) and the ex- 
treme value a = -1.73. We remind the reader that such extreme 
a value, that we think is not optimal for our da ta at this redshift, 
never t heless has been fo und in the literature dReddv & Steidell 
120091 ; lOesch et"ai] 1201 Ol) at very close redshift. Other findings 
from other works span at this z the range -1.7 < a < -1.1. 
In Table IB. fl we give a complete list of the LDs for our LFs de- 
scribed in Sec. 13.21 i.e. including the best LFs, those with all free 
STY parameters at z ~ 0.1 and z > 1.7, and with a set to -1.1 or 
-1.73 at z> 1.7. 

In the literature the treatment of cosmic variance as a source 
of uncertainty differs from work to work. If we consider the data 
sets overplotted in Figs.[3]and|5j we note the following. Some au- 
thors do not take into account the role of cosmic variance (e.£ 



ISchiminovich et al.ll2005l IWilson et al.ll2002l llwata et al.ll2007 
or co nsider it negligible with respect to statistical uncertainties 
(e.g JPaitani et al.l2007l) : o r consider it negligible be cause of the 
large volume explored (e.g. Ivan der Burg et al.ll20l"oL 4 deg 2 ); or 
mitigate its effects using several fields in different lines of sights, 
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Fig. 3. Left panel. VVDS Deep+Ultra-Deep FUV-band luminosity densities derived from our best LF determination (FUV-150 red 
filled circles) as a function of redshift. Error bars are explained in Sect. [4] Gray data points are ultraviolet LF estimates found in 
the literature, as detailed in the labels (s ee also at the end o f the caption). Right panel. FUV-band luminosity densities restricted 
to galaxies brighter than L* at z = in IWvder et all d2005l) (M FUV = -18.12). Symbols are like in the le ft panel. Black ope n 
circles and squares represent the LD down to Mpuv = -18.12 summing the 1 /V„ mx points in this work and in lTresse et al.l d2007l) . 
respectively, instead of integrating the Schechter LF fit. We did this only up t o z ~ 1 because in V VDS data at z > 1 the bias limit 
for l/Vmax computation is brighter than Mfuv = -18.12. List of references: iTresse et al] (120071) (VVDS Deep, FU V-150, empty 
upside -down triangles); compilation of FUV and NUV data extracted from iHopkinsI ( 2004 ) (thin open diamonds); ISteidel et al.l 
d 19991) (FUV-1700, fille d upside-down t riangles); ISchimin ovich et al] J20051) (F UV-150, filled squares th ose at z < 1.2 GALEX- 
VVD S forz < 1.2 data);lBouwens et all d2007l) (FUV- 1600, thick open squares): iReddv & Steidel d2009h (FUV-1700, filled trian- 
gles);|Wvder etal] J2005I) (GALEX-2dFGR S, FUV-150, th i ck ast erisk); |PaltMiLet_^ d2007|) ( FUV-1700, th i n aste risk) ;|l wata et al.l 
j2007D (FUV- 1700, thick open diamonds); iGabasch et al] d2004l) (FUV-150, open triangles); lOesch et al] d2010h (FUV-150, thin 
open squares): ISawicki & Thompson! d2006l) (FUV-1700, crosses). 



although gen erally these fields are much smaller than the VVDS- 
2h field (e.g.. ISteidel et~ai]ll999MReddv & Steidelll2009h : or the 
cosmic variance has been included as uncertainty in the LF nor- 
malisation (e.g jBouwens et aT]l2007l) . We choose to compute it 
with the most conservative method (see above) and to include it 
in our comprehensive uncertainty on the LD. 

The right panel of Fig. [3] shows the FUV-band LDs for 
galaxies brighter th an L* at z ~ 0, i.e., the determination by 
IWvder et alfd2005l) . corresponding to M* = -18.12. In Fig. [J 
we overplot data points of other studies which published the 
three Schechter parameters of their rest-frame ultraviolet LFs, 
that we integrated also to L/ ni „, = 10 15 W Hz -1 (left panel) and 
to L faint = U{z = 0) (right panel). 

From the left panel in Fig. [3] it is evident that the VVDS 
LD derived in this work shows a peak in its evolution with time. 
After an increase by a factor of ~ 2 from z ~ 4 to z ~ 2, the 
LD decreases sharply by a factor of ~ 4.5 down to z ~ 0. Fig. [2] 
shows that the LD is increasing from z ~ to z ~ 2 because 
of the continuous brightening of L* , with the normalisation and 
the faint end slope not changing so much. The brightening of L* 
also assures that the LD at z ~ 2 would be robustly larger than at 



z ~ 1.5, even if at z ~ 2 we had used a flat slop^l In contrast, at 
z > 2 the LD evolution is mainly driven by the evolution of the 
normalisation: the LD decreases from z ~ 2 to z ~ 4 because of 
the steep decrease of (p* . We note also that the decrease of the LD 
for z > 2 with our best estimate (red circles) is even smoothed 
away because we are assuming a increasing with redshift (that 
should cause the LD to increase). Only by using a much sharper 
steepening of a at z > 2 we could make the LD peak disappear, 
making the LD staying roughly constant at z > 2. For any other 
choice (the slope slowly increasing with z, or fixed to a flat or 
steep value like in the two extremes that we consider in the error 
bars), the LD is clearly decreasing at z > 2. 

Our result is an improv ement with respec t to the similar 
study of the FUV-band LD in lTresse etal] d2007l) , that was based 



on our first VVDS-Deep spectroscopic sample ( Le Fevr e_etal 



2005b) and our UBVRI photometric sample dLe Fevre et al 



l2004bl) . On one side, the so-called redshift desert (2 < z < 3) 
could not be analysed, as it is now thanks to our Ultra-Deep 
data. On the other side, the spectroscopic and photometric data 
were too shallow to determine a of the FUV-band LF, and thus 
a constant value (a = -1.6) was used over the explored redshift 



3 We remind that the lower limits of the error bars at z > 1.7 include 
the possibility of a = -1.1. 
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range (0 < z < 2, 2.70 < z < 5), as previously done in the 
literature (see Table IB. 2"b . Thanks to our new spectroscopic and 
photometric data, with the present work it is the first time that 
the presence of a peak at z ~ 2 in the FUV-band LD evolution 
has been robustly assessed from a single survey. 

4.2. Discussion and comparison with literature 

A large scatter between the different LD determinations is 
present in both panels of Fig. [3] In integrating a Schechter LF 
down to the faintest luminosities, the value of the faint-end slope 
has a strong impact, especially if it is steep, leading to very dif - 
ferent LDs for different a (see, e.g., Fig. 1 in lTresse et al.ll200"7l) . 
Table IB. 21 shows that the a determinations span a large range 
of values at every redshifts, and it is evident that our determi- 
nations of a are lower than the average, especially at z < 1.7. 
This is the main reason why our total LD appears on average 
lower than the other values in the literat ure (but still in agree - 
ment with works that used or — 1, like in lGabasch et al.ll2.Q04h . 
We have already discussed in Sect. [3] the difficulty to compare 
rest-frame FUV LFs, and thus LDs, from various surveys, con- 
trarily to rest-frame optical LFs and LDs. By integrating the LFs 
down to Mfuv = -18.12 (right panel of Fig.0, the differences 
in the values of a should have little impact on the resulting LD, 
when this integration limit is very close to M* (i.e. at z < 1). 
The observed scatter in the right panel of Fig. [3] means that L* 
and (p* are also discrepant among various data sets. Furthermore, 
since the three Schechter parameters (L*, <p*, a) are strongly cor- 
related, different values of (L*, 0*) arise even within the same 
data set when a is fixed to different values (see below). For in- 
stance, taking a steeper a leads to a fainter L* and a lower <t>* to 
fit exactly the same data set. 

In the right panel of Fig. [3] we also show our L > L* =0 
LDs obtained by summing the 1/V max values for galaxies with 
U < L < Li, r igh t (black empty circles). The discrepancies 
between this estimate and the one using the integration of the 
Schechter function mirror the level of extrapolation made using a 
STY fit. We can see that, within the error bars, our extrapolation 
is small, which ensures the goodness of our fit at Mfuy < -18.2 . 
The same exercise is repeated for the data in lTresse et al.l d2007l) 
(black empty squares). 

We can test the effect of different LF parameterisations on 
tota l LD estima t es, us ing two works within the VVDS sur- 
vey^ Tresse et all d2007l) used the VVDS Deep data presented in 
iLe Fevre et alJ d2005bl) (UBVRI) to derive the rest-frame FUV- 
150 LDs. Their slop e a has been s e t to — 1 .6, based on GALEX- 
VVDS results from lArnouts et alJ (120051) at 0.2 < z < 1, since 
the depth of the U-band data available at that time did not enable 
them to constrain the faint-end slope as we can do with our new 
u*g'r'i'zJHK s data. Therefore, their LDs at z ^ 1.5 are higher 
than those found in this work for both the total galaxy sample 
and the L > L* =Q galaxy sample. Nevertheless, if one compares 
the minimal LPs (i.e., t he sum of the 1 /V mflA ) between this work 
and lTresse et al.l (120071) . they are in fairly good agreement. 

Another example is given by the total LDs from 
iReddv & Steidell d2009t) at z ~ 2 and ~ 3, that are a factor of 
~ 2.5 larger than ours. The difference is mainly explained by the 
much steeper slope in their LF (see Table IB. 21 . In fact, from the 
left panel of Fig. [3] we see that their LD is in much better agree- 
ment with the upper limits of our error bars, that correspond to 
the extreme case with a = -1.73. 

It is worth noticing that lOesch et alJ (1201 Oh find a constant a 
within 0.5 < z < 2.5, but much steeper than ours (a ~ -1.7). 
Moreover, they do find a peak in the LD evolution as we do, but 



at z ~ 2.5. We remark that they have to use points from the liter- 
ature to find this peak (the highest z of their LD determinations 
being z = 2.5), while we are able to constrain the LD peak using 
our own homogeneous data. 

More generally, comparison with literature appears tricky, 
not only because of very different assumptions on the LF shape 
(the need to set a or M*, for instance), but also because the LF 
shape can be intrinsically different according to the various sur- 
vey characteristics (flux limit, band of selection...). Due to this 
difficulty in comparing different works, it is clear the importance 
of tracing the rest-frame FUV LDs over a wide range of redshifts 
with a single and homogeneous survey. This is what we have 
been able to do with our data set. 



5. The dust attenuation and the star formation rate 
densities up to z - 4.5 

A robust determination of the dust attenuation is a critical 
element to transform the luminosity density into the effec- 
tive star formation rate density. Thanks to our present VVDS 
Deep+Ultra-Deep spectroscopic sample, combined with a wide 
range of deep multi-wavelength data (u*g'r'i'zJHK s ), we can 
now study in details and with homogeneity the cosmic dust at- 
tenuation evolution over < z < 4.5. 

In the following subsections, we explain how we determined 
the dust attenuation in our sample and how it compares with (and 
improves upon) other measurements in the literature (Sec. l5.lt . 
Then, we show how we derived the dust-corrected SFRD for 
the global galaxy population (Sec. 15.21 1. and how galaxies with 
different luminosities contribute to the total SFRD CSec. l5.31 >. 

5. 1 . The dust attenuation 

Estimating the average dust attenuation properties in a galaxy 
survey is a complex task. A number of authors have suggested 
various methods to derive the dust attenuation using a multi- 
wavelength data set. This leads to one important source of scat- 
ter in the literature in the determination of the dust corrected 
SFRDs. For instance, dust attenuation can be estimated using 
dust reddening curves or UV slopes in the FUV, or the Balmer 
decrement measured from Ha and H/3 (with the drawback that 
Ha is difficult to observe for large samples at z > 0.3). It can be 
also derived with multi-wavelength SED fitting techniques, es- 
pecially when the knowledge of spectroscopic redshifts enables 
to break a dust-redshift degeneracy. Nevertheless, if there is a 
large and prominent population of highly dust-obscured galax- 
ies, it will be missed by deep rest-frame FUV data sets. In prin- 
ciple, a powerful method is based on adding the total SFR ra- 
diation, that is the one passing through dust (FUV emission) 
and the one re-emitted by dust (FIR emission). However, this 
method has potential drawbacks. For example, results will be 
exclusively linked to those galaxies visible at both wavelengths, 
and mainly at z < 2. Furthermore, the situation in infrared sur- 
vey observations is not optimal due to unknown contents of cold 
and hot dust, or broad-line active galaxy nuclei contamination, 
or large PSFs making difficult the measurement of uncontam- 
inated s i ngle s ources. We refer the reader to the discussion in 
ICalzettil J2008h for details on computing dust attenuation using 
both FIR and FUV data. 

Since we want to derive the dust attenuation in a similar man- 
ner over < z < 4.5, we estimate the dust attenuation using the 
results from the SED fitting analysis. This method is not always 
fully reliable on a object-by-object basis. Still, it provides a good 
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Fig. 4. Dust attenuation Afjjy in magnitudes as a function of redshift. The thick solid red curve represent the fit (see text for details) 
of the average dust attenuation, as a function of redshift, determined in this work, using the Calzetti's law. The Afuv values derived 
from this fit in eac h redshift bin are lis ted in Table Q] The magenta dot-dashed curve is the average Afuv determined in this work 
using the recipe in lMeurer et al.l (1 19991). b ased on the B slope. The green dotted curve is the same as the dot-dashed one, but here 
we use the recipe in ICortese et al.l (120061) . calibrated with normal star forming galaxies. The thin blue horizontal solid line is a 
constant Afuv computed with the Calzetti's law using one single typical value of E(B - V) (=0.13). The right y - axis shows 
the multiplicative factor to be applied to the observed luminosities, i .e., lQ a44f »' / All th e other symbols represent dust attenuations 
found in literature, as indicated in the labels: orange dashed line from lTresse et alJ d2007l) (obtained comparing not-corrected SFRDs 
with dus t -corre cted SFRDs from other works, see text for details); gray open diamonds and filled squares from the compilation in 
iHopkinsI (12004 H04) (only FUV and NUV determinations, diamonds and square s for a correction independent from or dependent 
on SFR, respectively); cyan filled upside- down triangles fr om Steidel et al.l (Il999[ St99): thick gr ay crosses from lOuchi etal] (12004 
O04); orange open and filled squares from lBouwens et aL (2009, B09) and lBouwens et al.l d201 ll Bl 1) respectively (s lightly shifted 
in red shift for clarity); blue filled circles from lTakeuchi et all d2005l T Q5). based on the 02h field GALEX-VVDS LFs (lArnouts et all 
120051) and the CDFS field Spitzer 15-yum LFs dLe Floc'h et al.ll2005T) at 0.2 < z < 1; arr ows fromlSchiminovich et all (120051 Sc05t. 
also based on GALEX-VVDS survey; filled magenta triangles and magenta crosses from lReddv & Steidell d2009l R&S09) (triangles 
and crosses for a dust attenuation dependent on or independent from FUV luminosity, respectively). 



estimate of the average attenuation properties of the whole popu- 
lation under study and of their evolution with time. Also it avoids 
the use of restricted cross-matched optical-FIR catalogues. 

We primarily use the recipe in ICalzetti et alJ d2000t) for ac- 
tively star-forming galaxies to derive our dust attenuations at 
1500 A, and we also report and discus s those derived with other 
recipes. We follow the prescription in lCalzetti et alJ d2000l) (see 
Eq. 4), that is, A{A) = E(B - V) star k(A), where E(B- V) srnr is the 
intrinsic colour excess of the stellar co ntinuum of a galaxy, and 
k{A) is the starburst reddening curve in lCalzetti et alJ d2000 fl at 
[0.12-0.63]//m. As described in Sect. 12. 21 our template SED fit- 
ting assigns to each galaxy a value of E(B - V) srnr chosen from a 
grid of five possible values (0.1,0.2,0.3,0.4,0.5). In each stud- 
ied redshift bin, we compute our resulting mean E(B - V) star and 
A fuv ■ We remark that our E(B-V) s t ar computations do not seem 
to depend on observed (i.e. not dust-corrected) rest-frame FUV 
luminosities, which indicates that our sample includes a large 
variety of actively star-forming galaxies at each luminosity. For 
this reason, we are confident that our mean E(B - V) slar values, 
obtained averaging galaxies brighter than the given magnitude 
bias at each redshift, do not depend on the range of luminosity 
spanned at each redshift. 



k(X) = 2.659 x (-2.156 + 1.509//1 - 0.198//1 2 + 0.01 l/A 3 ) + 4.03 



Fig.|4]shows our results. We find that Afuv increases steeply 
and fast from z ~ 4 to z ~ 2, i.e. of 1 mag within only 2 Gyr. It 
keeps on mildly increasing from z ~ 2 to z ~ 1 (~ 0.2 mag in 
2.5 Gyr) and then it continuously decreases by ~1.1 mag within 
~ 8 Gyr from z ~ 1 to z ~ 0. This results in a maximum dust at- 
tenuation value at z ~ 1 (Apuy 2.2). We smooth Afuv{z) with 
a linear fit and a 3 rd order polynomial function for the redshift 
ranges z < 1 and z > 1, respectively, to preserve the different 
trends of Afuviz) below and above this redshift. The red solid 
line in Fig. [4] shows our smoothed Afuy, and we give its values 
in each redshift bin in Table Q] We remark that the chosen fits 
depart from the fitted data by a very small amount, so we do not 
overplot the original Afuy for clarity. 

We also comput e our average Afuy with the recipe given by 
iMeurer et al. based also on actively star-forming galax- 

ies, but using the B slope of the UV continuum as a proxy to the 
UV attenuati on (dot-dashed mag enta line in Fig. [4]). We use the B 
expression of lKong et al.l d2004l) which accounts for the star for- 
mation history to reduce the large scatter of this method. Using 
the B slope l eads t o larger Afuv, in particular at 1 < z < 2. 
iTrever et al.l d2007l) show that this method works well at least 
up to z = 1, but breaks down for red-sequence and very blue 
compact galaxies and for the majority of high-z Lyman break 
galax ies, which form a low-a ttenuation sequence of their own. 
Also, IWiiesinghe et alJ d201 ll) emphasise that this method over- 
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estimates SFRs unless a modified relation between ft and the 
attenuation factor is used as indicated by other galactic prop- 
erties than the UV continuum, in particular for high-z galaxy 
SFRs. Therefore, a good knowledge of the galaxy population un- 
der study is necessary for a reliable application of this method. 
For instance, we comp ute our average valu es of Apuv with the 
/? expression given by ICortese et al.l (120061) . calibrated on nor- 
mal star-forming galaxies, and we plotted them in Fig.|4]as well 
(green dotted line). While this relation follows the same trend 
as the previous estimates, the resulting attenuations are lower by 
more than 1 m ag. On the ba s is of t his analysis, we conclude 
that taking the ICalzetti et al.l d2000l) reddening curve appears 
the best solution for the VVDS Deep+Ultra-Deep sample since 
it is largely dominated by actively star-forming galaxies (from 
~ 40% at z ~ 0.1 up to ~ 85% at z ~ 4.0, see Appendix [B]). 
Nevertheless, whatever the method to compute the dust attenua- 
tion, the amount of dust increases from z = 4 to z — 1, and then 
decreases down to z ~ 0. 

In Fig. |4] we also overplot Apuv dust attenuation values 
taken from the literature. We can see that data are diversely 
spread between 1 and 1.8 mag. Actually, they are often de- 
rived in d ifferent ways, dependin g sometimes on FUV-band lu - 
minosity dReddv & Steidelll2009i). or on SFR dHopkinsI 120041) . 
or being a constant value dSc himinovi ch et al. 112005] GALEX- 
VVDS). Conversely. iTresse et al.l d2007b derived them by c om- 
pari ng the dust corrected 12^m dPerez-Gonzalez et al.l2005l) and 
Ha dTresse & Maddox|[l998T: ITresse et al.ll2002l) SFRDs to their 
uncorrected FUV-derived VVDS Deep SFRDs. They found that 
the attenuation at 1500 A was ~2 mag from z — 2 to z — 0.4, 
and then it decreased from z = 0.4 to z — down to ~1 mag. 
These last Ap yy values are very close to our present work. We 
plot them in Fig. [4] with an orange dashed line. We are in good 
agreement also with iTakeuchi et al.l d2005l) at 0.2 < z < 1 . To 
estimate the mean dust attenuation , they compared the GALEX- 
VVDS LFs dArnouts et al.l 120051) to the Ch andra Deep Field 
South Spitzer 15-/U LFs dLe Floc'h et al.l2005l). O ur Afhv agrees 
also with the one deriv ed by Steidel et alJ ( 1999|) at z ~ 3 and by 
iBouwens et al.l d2009l) at z ~ 4 Tbut see lBouwens et alJl201 ll for 
a much lower Apuy value). 

In summary, with the comprehensive VVDS Deep+Ultra- 
Deep sample we are able to compute the evolution of the mean 
cosmic dust attenuation over ~12 Gyr. It increases rapidly of 1 
mag in 2 Gyr (from z ~ 4 to z ~ 2), then it reaches a maximum 
at z ~ 1 (increasing by ~ 0.2 mag in 2.5 Gyr), and finally it de- 
creases continuously down toz ~ 0(~1.1 mag within ~ 8 Gyr.) 
There is a maximum in dust attenuation at z ~ 1 (Apuy - 2.2). 



5.2. The dust-corrected SFRD 

To transform FU V fluxes into star for mation rates, we use the 
SFR calibration of lMadau et all dl998l) . It yields: 



SFRD(z) = 1.4 10" 28 LD FUV (z) 10' 



0.4 A FUV (z) 



(1) 



where the SFRD is in M Q yr _1 Mpc _ units and the LD in erg 
s _1 Hz 1 Mpc~ 3 . This formula i n clude s the dust attenuation 
Afuv(z), and assumes a ISalpetetl d 19551) initial mass function 
(IMF) including stars from 0.1 to 125 solar masses. The resulting 
dust-corrected SFRDs are shown in Fig. [3] The values of uncor- 
rected (i.e., Afuy{z) = 0) and dust-corrected SFRDs are given 
in Table Q] For the SFRDs, we assume the same uncertainties as 
the LDs (see Sect.O, that is, we do not include any dust redden- 
ing errors since it is the choice of the method which dominates 
(see Sect. 15.11 ). Still, we include STY fit errors, Poisson noise, 



cosmic variance and weighting scheme errors. We remind that 
the STY fit errors, when a is fixed, include the span in LD (and 
so in SFRD) that one would have setting a at the two extreme 
values (—1.1 and -1.73) discussed in Sec. [3] We overplot in the 
top panel of Fig.[5]dust-correc ted UV-derived SFRDs from other 
published works, assuming a ISalpeterl d 1955b IMF. We refer to 
iHopkins & Beacoml(l2006l) for a detailed comparison among dif- 
ferent IMFs. 

Like for the LD, we find a peak in the SFRD evolution at 
z ~ 2. This peak is preceded at earlier cosmic epochs by a rapid 
increase of a factor 6 from z ~ 4.5, then followed by a decrease 
by a factor of 12 to z ~ 0. We note that the VVDS Deep+Ultra- 
Deep SFRDs are globally, within error bars, in agreement with 
the literature of UV data sets obtained at various redshifts, but 
it is the first time that the cosmic SFRD history is continuously 
traced from z = 0.05 to z = 4.5 in a homogeneous way within 
one sample. We can therefore assess the presence of a peak in the 
SFRD evolution at z ~ 2, without speculating about the different 
selection functions of disparate surveyfl We fitted the SFRD as 
a function of z below and above z = 2, using SFRDoc (1 + zf. 
We find/? = 2.6±0.4atz < 2, and/? = -3.6 + 1.9atz > 2. These 
two fits are shown in the bottom panel of Fig. [5] 

We remark that this peak is the result of the combination of 
the LD evolution and the dust attenuation evolution, but it is not 
uniquely driven by the dust correction. First of all, we already 
found a clear peak at z ~ 2 in the evolution of the LD. Secondly, 
the peak of the dust attenuation evolution does not coincide, in 
terms of time, with the SFRD peak, because we find it at z ~ 
1 and not at z ~ 2. Surely, the shape of the dust attenuation 
evolution (steeply decreasing from z ~ 2 to z ~ 4) enhances the 
relative decrease of the SFRD with respect to the decrease of the 
LD in the same redshift range, but it is not creating the SFRD 
peak by itself. 

Although the scatter of measurements beyond z ~ 1 is in- 
creasingly larger, at z > 3 our measurements are lower than most 
other measurements. This is mainly due to our use, at z ~ 3, 
of a flatter faint-end slope in the LF than other literature stud- 
ies, and to a very low <p* at z ~ 4 (see Sec. 13.3.21 for a discus- 
sion about <p* at z > 2). The slope of the LF remains a major 
source of uncertainty at these redshift s, but also the computa tion 
of the dust attenuation. For example, IBouwens et al ]d201 ll) re- 
cently revised their previous Apuv computation dBouwens et all 
120091) : at z ~ 4, their dust attenuation correction factor (com- 
puted for L > 0) is now ~ 40% lower than before. This leads 
to a lower total dust-corrected SFRD by ~ 0.2 dex, giving fur- 
ther sir2£ortJo_ojrrre^uIt of a low SFRD at this redshift (see 
also lCas t ellano et aUl201 ll). In terestingly, the dust attenuation at 
z ~ 4 in IBouwens et al.l d2009l) was closer to our Ap ay, but their 
new determination of SFRD is closer to our SFRD value. 

In the bottom panel of Fig. [5] we add also SFRD data points 
in the literature, derived using other SFR calibrators than the UV 
one, mainly at z < 2. The non FUV-derived SFRD estimates do 
not increase the scatter observed with the FUV-derived ones. At 
z < 2, the cosmic SFRD evolution can be also derived from the 
stellar mass density assembly, that is from the integrated history 
of the SF RD down to a give n epoch. For instance, we take re- 
sults from lllbert et al l ( 1201 Ob (their Tables 2 and 3), who have 
computed the stellar mass function s and p stn r densities o ver the 
multi-wavelength COSMOS field dScoville et al.l 120071) . From 



5 Note that also uncertainties are often computed in different ways 
in each galaxy survey, and they are sometimes underestimated because 
not all-inclusive of many sources of errors. This further increases the 
non-homogeneity of the SFRDs determinations. 



11 



O. Cucciati, L. Tresse, O. Ilbert et al.: The CSFRD since 12 Gyr traced by the VVDS 



z = 0.2 to z = 2, in each redshift bin we obtain p star densi- 
ties in summing the individual stellar mass densities determined 
for the "quiescent, "intermediate activity" and "high activity" 
galaxy populations. Assuming a linear fit, we recover its evo- 
lution with redshi ft, Pstar(z), and the im plied SFRD(z) with the 
recipe detailed in IWilkins et all d2008lrl Note that we assume 
that 30% of the created stella r mass re t urned into the interstel- 
lar medium: as suggested by iPrantzosI d2008l) . this is the aver- 
age fraction derived using different IMF. As shown in Fig. [5] 
the resulting inferred SFRD(z) agrees on average with our di- 
rect SFRD(z) estimate. Nevertheless, there is some discrepancy 
(~ 0.2 dex) around z ~ 1, which happens at the cosmic dust 
attenuation peak (see Fig.|4]i within the star-forming population. 

In summary, the VVDS Deep+Ultra-Deep survey traces con- 
tinuously the cosmic dust-corrected SFRD over ~12 Gyr in 
a comprehensive way, with homogeneous treatment of data, 
sources of errors, and dust attenuation correction. A clear peak 
of the SFRD emerges at z ~ 2, preceded at earlier cosmic times 
by a rapid increase by a factor ~6 within ~2 Gyr, and followed 
by a decrease by a factor ~12 within ~10 Gyr, down to z ~ 0. 
For the first time, the SFRD evolution is genuinely established 
over 12 Gyr. 

5.3. The contribution of galaxies with different luminosities 

The relative contributions to the total SFRD of distinct galaxy 
populations vary with redshift, an d it is their combin ation that 
shapes the global SFRD evolution. iTresse e t al. (2007]) have as- 
sessed the contributions of galaxies with different luminosities 
to the non dust-corrected FUV-band LDs, showing, for instance, 
that at z ~ 3.5, the different FUV populations equally contribute 
to the total LD. Then, as the Universe is ageing, the contribution 
of the luminous FUV population {Mfuv < -19) dimes, while 
the one of the fainter population increases. With our present data, 
we can directly study the contribution of galaxies with different 
luminosities to the dust-corrected SFRD. We defer to a future 
work the analysis on how galaxies contribute to the total SFRD 
as a function of their stellar mass. 

Fig. [6] shows the dust-corrected SFRD of galax- 
ies with different luminosity upper limits: Mfuv ^ 
-17.5,-18.5,-19.5,-20.5,-21.5 and M FUV < M* FUV (we 
remind that M* FUV varies with z, see Table[T]i. We apply to each 
SFRD the same dust correction derived for the total SFRD, since 
we do not detect any significant E(B - V) colour dependence 
with the FUV luminosity (see Sect. 15. 11 1. We see that, as times 
goes by, the total SFRD is dominated by fainter and fainter 
galaxies. In fact, bright galaxies with Mfuv ^ -21.5 contribute 
little to the total SFRD for z < 2, and this holds for galaxies 
with M FUV < -20.5 at z < 1 and with M FUV < -19.5 at z < 0.8. 

To investigate further these trends, we made 6 classes of 
galaxies within the following Mfuv luminosity ranges: > -17.5, 
[-17.5,-18.5], [-18.5,-19.5], [-19.5,-20.5], [-20.5,-21.5] 
and < -21.5 (see Fig. |7J. We describe below the three main 
cosmic eras of the SFRD history. 

From z = 4 to z = 2. During this period, the SFRD from 
each galaxy subclass increases, with the only exception of the 
brightest galaxies {Mpuv < -21.5), whose SFRD starts decreas- 
ing at z ~ 3. Even including all sources of uncertainties, the 
contribution of these brightest galaxies to the total SFRD is a 
maximum of 20% at z ~ 3 and decreases to a maximum of 10% 
at z = 2. The net effect is that the total SFRD increases from 

6 St ellar masses in lllbert et al .1 j2010h are comp uted using a | Chabriej 
d2003h Initial Mass Function, that we convert to a lSalpeteJdl955h IMF. 



z — 4 to z — 2, because the SFRDs of the main contributors 
{Mfuv > -21.5) increase as well. 

From z-1toz- 1 . The period corresponds to the change of 
the shape and normalisation of the FUV LFs (see Figs. [1] and [SJ. 
Thus, the trends are more differentiated from class to class: the 
SFRDs for [-19.5,-17.5] stay about constant, for > -17.5 and 
[-20.5,-19.5] mildly decrease, for < -20.5 sharply drop. The 
net effect is however a global decrease, the highest star-forming 
populations being shut down (galaxies brighter than Mfuv = 
-20.5 contribute maximum by 3%). 

From z = 1 to z ~ 0. Through this period, the individual 
SFRDs reach a general decreasing trend, particularly steep for 
[-20.5, -19.5] galaxies. As times goes by, the less galaxies are 
forming stars, the larger is their contribution to the total SFRD. 
The net effect is that the total SFRD decreases from z = 1 to 
z~0. 

To summarise, the total SFRD evolution (with an increase 
from z — 4 to z — 2, a peak at z — 2 and a decrease from z = 2 
to z = 0) is driven by different galaxy classes at different epochs. 
In particular, the percentage contribution to the total SFRD of all 
galaxy classes is roughly the same at z = 2, but the peak of the 
total SFRD at this redshift is mainly due to a similar peak of the 
SFRD from galaxies with -21.5 < M FUV < -19.5 (L > L* =2 ). 
In contrast, the contribution of the most luminous star-forming 
galaxies reaches its maximum at higher redshift (z ~ 3). We note 
that the faintest galaxies {Mfuv > -17.5) are the only ones that 
show a continuously decreasing SFRD from z = 2 to z = 0. 

Fig.Qshows that at z > 3 each FUV luminosity class shows 
a very similar SFRD trend with z, while for z < 3 the contri- 
butions of the various classes to the SFRD become more and 
more dissimilar as times goes by. This implies i) the presence 
of a very high-z population not yet strongly affected by star for- 
mation regulation/quenching mechanisms and ii) a downsizing 
pattern in the SFR history. 

We remark that the contribution of the faintest galaxies 
{Mfuv > -17.5) is clearly affected by the shape (and its evo- 
lution) of the LF. We refer the reader to ITresse etal] d2007l) for 
an extensive analysis of the role of a in the relative contribution 
of faint galaxies to the total LD. The difficulty in constraining 
a, and the spread in typical luminosity of the targeted popu- 
lations, are surely two important sources of uncertainty in the 
literature in the measurement of the total SFRD shape in this 
redshift range. 



6. Summary & Discussion 

In this work, we have computed the rest-frame FUV luminosity 
functions, luminosity densities and dust-corrected star formation 
rate densities within a single and deep redshift survey, merging 
the VVDS Deep and VVDS Ultra-Deep data sets (with over- 
all 17.5 < Iab ^ 24.75), over a 12 Gyr cosmic time baseline 
(0.05 < z < 4.5). We also derived the average dust attenuation in 
the entire redshift range explored. Our data constitute a large im- 
provement with respect to those used in previous VVDS studies: 
i) the Deep survey includes more redshifts, ii) repeated observa- 
tions (leading to 100% secure redshift measurement for previ- 
ous lower quality spectra) enabled us to correct our «(z) for flags 
1 and 2, iii) we use the Ultra-Deep sample, pushing 0.75 mag 
deeper, iv) we have deeper optical photometry, in particular in 
the M*-band and z-band, and v) we have used new near-infrared 
deeper photometry. 

Our results are summarised as follows. 
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Fig. 5. Top panel. Total dust-corrected UV-derived SFRDs as a function of redshift from the VVDS Deep+Ultra-Deep sample (red 
filled circles, see Tabl e [p. Our uncertai nties are explained in Sect. [3] The black dashed line is the SFRD(z) implied from the 
stellar mass density in lllbert et al.l (l2010h (see text). We overplot other results from the literature, as detailed in the labels and at 
the end of the caption. Our SF RDs, as well as thos e from the literature, are derived using the FUV-band L Ps converted in to SFRD 
with the scaling relation from iMadau et al.l (Il998l) . All data have been homogenised with the same IMF (I S alpeterl 1 1 95 5h . Bottom 
panel. Total dust-corrected UV-derived SFRDs as a function of redshift from the VVDS Deep+Ultra-Deep sample (red shaded area, 
corresponding to red circles and error bars in the top panel). The two long-dashed lines represent two fits to our SFRDs in the form 
oc (1 +zf (see text for details). Gray points are from literature, derived from various SFR calibrators (UV, emission lines, IR, radio): 
th ey includ e the li terature points as in the top panel, plus other works as in the labels. L ist of references, top panel. Compilation 
in iHopkinsI (|2004|) taking on ly the FUV and NU V determinations (grey open diamonds); ISteidel et al.l (119991) (LBG sa mple, cyan 
solid upside-down triang les): lOuchi et al.l (12004 (S DF and SXDF LBGs sample, thick grev crosses'): I Wvder et all (120051) (GALEX- 
2dFGRS, bold asterisk ); ISchiminovich et al l d2005h (GALEX-VVDS at z < 1.2 and HDF above, small so lid squares with a r rows) ; 



iBouwens et al. | (120091) and IBouwens et al.l d201 llT ( LBG s, orange open and filled squar es, respectively): iReddy & Steidell d2009l) 
(LBG, solid magenta triangles: Ivan der Burg et all (1201 Ol) (CFHTLS, black open circles'): ICaste llano et al.l d201 ll) (LBG. blue open 



circles). Li st of references, bottom panel. All the po i nts in the top panel, plus: grey thin diamonds comes from the entire compilation 
report ed in Hopkins ( 20041); Perez-Gonzalez et al.l d2005l) (12/itn, upside-down and normal solid triangle s are the lower and upper 
limits: iRodighiero et al.l d2010l) (LI R Spitzer VVDS-SWIRE & GOODS, solid stars'): iBardelli et al.l d2010l) (radio VLA-zCOSMOS, 
open stars); iLe Floc'h et al. ni2005l) (24//m Spitzer CDFS, solid line; from their Fig. 14). 



- We find a flat and constant faint-end slope in the FUV-band 
LF at z < 1.7 (a ~ 1). We verified that this is unlikely to be 
the result of missing faint galaxies from our /-band selection, 
and that dust may have a role (see discussion in Sec. I3.3t . At 
the same time, M* FUV increases by ~ 1 .5 mag from z ~ to 
z ~ 1.2, while <p* starts decreasing at z > 0.7 (more than a 
factor of 2 decrease from z ~ 0.7 to z ~ 1 -5). 



- At z > 1 .7, we set a evolving with (1 + z), and this way it be- 
comes as steep as -1.73 at z — 4, consistent with values from 
deep photometric studies. In the meanwhile, M* FUV keeps on 
brightening (by an added 2.5 mag up to z ~ 4), and (p* de- 
creasing (by another factor of ~ 40 up to z ~ 4). We find 
that at z — 2, 3, 4 our M* FUV is on average brighter than what 
has been found in previous works, and that <p* is on average 
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Fig. 6. The dust-corrected SFRD from galaxies with different 
luminosity limits, as described in the labels. The background 
gray shaded area is the total SFRD, while the SFRD for dif- 
ferent luminosity cuts decreases for brighter galaxies. The lumi- 
nosity limits are kept fixed at all redshifts. The shaded area for 
each SFRD accounts for the same kind of errors as for the total 
SFRD (STY ellipses, cosmic variance, Poisson noise, weighting 
scheme, and at z > 1.7 also the span of a between a = -1.1 
and a = -1.73). When the luminosity limit is fainter than the 
LF magnitude bias (see Sec. 13.11 ) at a given redshift, the area be- 
tween the lower and upper 1-cr error is not filled but left empty. 
The dashed black line is the SFRD from galaxies brighter than 
M* at each redshift (i.e., the luminosity cut varies with redshift 
following M*). At each redshift, error bars on this SFRD are 
similar to those for the cut in luminosity closest to M*, and we 
do not overplot them for clarity. We apply to all the points the 
same dust correction as for the total SFRD, as we do not detect a 
clear dependence of E(B-V) on FUV luminosity (see discussion 
in Sec. 15711 



smaller, in particular at z ~ 4. A summary raw scheme of 
these trends within < z < 4 is shown in Fig. [8] 

- We find that at z ~ 3, while the projected number counts of 
I-selected galaxies like VVDS are at least twice larger than 
the projected number counts of LBG-selected galaxies (see 
Le Fevre et al., in prep.), the VVDS LF has about 50% higher 
density at Mfuv ~ -22 than the LF derived from LBG 
counts. To interpret this apparent discrepancy, we verified 
the importance of the redshift distribution n(z) shape when 
transforming number counts (as a function of the observed 
magnitude) into a LF, or vice-versa. Also, when comparing 
number counts of two different samples in a given redshift 
range, one has to take into account the n(z) of the two data 
sets. 

- We derived the evolution of the dust-attenuation in the FUV- 
band (Afuv) in the range 0.05 < z < 4.5, using a SED fitting 
method, in a consistent way from a single survey with a well 
controlled selection function. We find a continuous increase 
of Afuv by ~ 1 mag from z ~ 4.5 to z ~ 1 (with the increase 
being very slow at 1.5 > z > 1), then a decrease by the same 
amount from z ~ 1 to z ~ 0. This is the first time that the 
Afuv evolution has been assessed homogeneously on such a 
broad redshift range. 



Look-bock time (Gyr) 
0.0 7.7 10.2 11.4 12.0 12.3 




z 



Fig. 7. As in Fig. [6] but this time galaxies are divided in nar- 
row luminosity bins, as indicated in the labels. The first and 
last bins are actually the extremes of the luminosity distribution 
(the faintest and the brightest galaxies). Also in this case, as in 
Fig.|6l the errors include all sources of uncertainty. The two lines 
for each luminosity bin define the region of +l-cr error. Not to 
crowd the plot, as the plotted quantities have larger error bars 
than in Fig. [6j here we do not distinguish when the luminosity 
bin considered is fainter or brighter than the LF bias. This can 
be inferred by Fig. [6] 



- We traced the dust-corrected SFRD evolution over the past 
~ 1 2 Gyr. It can be schematically fitted as SFRD(z) oc ( 1 +z~f, 
with /3 = 2.6 + 0.4 at z < 2 and/3 = -3.6 + 1.9 at z > 2 (see 
Fig.|5]and[8]l. Thanks to the homogeneity of our data over 
such a cosmic time, we have been able to unveil the presence 
of a peak at z ~ 2 in the cosmic SFR history. This peak is 
preceded by a rapid increase by a factor 6 from z ~ 4.5, 
then followed by a general decrease by a factor 12 to z ~ 0. 
We remark that the epoch of the peak of the dust-corrected 
SFRD (z ~ 2) does not coincide with that of the maximum 
of the dust attenuation evolution (z ~ 1), and that a peak at 
z ~ 2 is already present in the evolution of our LD. 

- Studying the contribution to the total SFRD of galaxy pop- 
ulations with different properties, we find that as times goes 
by, the total SFRD is dominated by fainter and fainter galax- 
ies. Moreover, the presence of a SFRD peak at z ~ 2 is 
mainly due to a similar peak within the population of galax- 
ies with -21.5 < Mfuv ^ -19.5, while the most ex- 
treme star-forming galaxies reaches their maximum activity 
at higher redshift. 

Our data therefore consistently show a peak in the LD and 
SFRD at z ~ 2, where the LF is well constrained. While the 
increase in SFRD from z — 4 to z = 2 is not in question, the 
exact amplitude of this increase remains to be investigated, as 
the faint end slope of the LF is still unconstrained even from 
our very deep spectroscopic survey. Other measurements of the 
SFRD beyond z — 2, using e.g. photom etric redshifts derived 
from multi-band imaging (for example [G abasch et all 120041: 
ISawicki & Thompson! l2006t iBouwens et al.1 120091) . go deeper 
than spectroscopic samples but carry significant associated un- 
certainties either on the photometric redshifts or because of the 
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Fig. 8. Scheme of the total dust-corrected SFRD as a func- 
tion of redshift, as found in this work (solid broken line, see 
Sec. I5.2I ). We distinguish four epochs, separated with vertical 
dotted lines. At each epoch, arrows indicate the behaviour of the 
three Schechter parameters and Af yy: arrows pointing up mean 
that the parameter is increasing with increasing redshift, arrows 
pointing down mean that the parameter is decreasing. Red ar- 
rows indicate trends that in principle should decrease the total 
SFRD, the opposite is for blue arrows. Thick arrows indicate 
which parameter is driving the SFRD evolution at that given 
epoch. In this raw scheme, we neither consider Afuv dominat- 
ing, nor making increase or decrease the SFRD: clearly, its trend 
enhances or weakens the behaviour of the LD when it is trans- 
formed to dust-corrected SFRD, but the peak at z ~ 2 is already 
found in the LD. 



small size of the fields and/or samples. It is then fair to say that 
the exact increase of the SFRD from higher redshift to z ~ 2 
remains to be robustly quantified. 

The correct determination of the shape of the SFRD evolu- 
tion is necessary to understand which physical processes mostly 
affect galaxy evolution. The SFRD is the result of the transfor- 
mation of gas into stars and therefore requires a significant gas 
reservoir to sustain a strong star formation rate for a long time. 
A number of processes are expected to modify the gas reser- 
voir hence the SFR, including the efficiency of star formation, 
cold accretion along the cosmic web filaments, mergers with 
gas-rich galaxies, stellar feedback, SN feedback blowing gas out 
from the galaxy core, AGN feedback, cosmic photoionizing ra- 
diation, or environment effects which ma y result in star f orma- 
tion quenching (e.g., among many others , White & FrenkJ 
Efstathiod fl99 2i: ICole et al.lll994l l2000t iDi Matteo et al 



1991 



2005 



Baughl l2006t ICox et al.ll2008t IDekel et al.ll2009t Ide Ravel et al.l 
20091) . The exact balance of these different processes along cos- 
mic time will result in the observed SFRD. The SFRD peak that 
we find at z - 2 is produced by galaxies with L close to or even 
brighter than the L* of that redshift, requiring that significant gas 
reservoirs still exist at this epoch and are probably replenished 
by cold accretion and wet mergers, while feedback or quenching 
processes are not yet strong enough to lower the SF. The knee 
of the rest-frame FUV LF shape is smoothed away at z > 2, i.e. 
there is not anymore a clear distinction between the high and in- 
termediate star-forming galaxy densities. It does not mean that 
the high-z gas-rich galaxies form a homogeneous population, but 
likely that the long time-scale (t > 3Gyr) processes do not heav- 
ily affect the galaxy SFR yet. Knowing the SFRD, we may hope 



to identify the relative contribution of these different processes 
at different epochs. 

Using simulations (semi-analytical galaxy evolution mod- 
els, smoothed particle hydrodynamics simulation), several 
authors have attempted to reproduce the observed Cosmic 
Star Formation History with theoretical predictions (e.g . 
[ Baugh et alJl2005b ISomerville et al.ll2008t iFontanot et al.l 120091: 
iHookins et al.ll2010Hvan de Voort et aLll201 ll) As a recent exam- 
ple. IWeinmann et all ( 201 1 ) use semi-analytical galaxy evolu- 
tion models to predict the SFRD evolution (see their Fig. 10). 
They start from a standard model of galaxy evolution, and pro- 
duce slightly different versions of it by tuning one or more ingre- 
dients at a time (such as star formation efficiency, stellar feed- 
back, merger processes...). Although the main goal of their work 
is not to reproduce the SFRD evolution, from their Fig. 10 it is 
clear that models in which different physical processes have been 
enhanced or depressed predict significantly different Cosmic 
Star Formation Histories. All of their models reproduce qual- 
itatively the fast SFRD increase from z ~ to z ~ 1, but at 
z > 1 the predicted SFRDs have different behaviours, showing 
growths and decreases at different rates, in some cases a plateau, 
in some cases a more or less pronounced peak at different possi- 
ble redshifts. It is worth noticing that one of the models (where 
feedback has been tuned) predicts a sharp peak in the SFRD evo- 
lution at z ~ 1.5, qualitatively similar to the one that we find at 
Z ~ 2. It is not the aim of this paper to compare in details our 
findings with model predictions, and we will address this issue 
in more details in a future work. Here we want only to stress 
the importance of bringing strong observational constraints on 
the SFRD from a unique and homogeneous galaxy sample cov- 
ering a large cosmic time of ~ 12 Gyr, which will need to be 
reproduced by next generation models. 

Another remarkable finding of our work is that the peak of 
the dust and the SFR do not coincide, differently from what one 
could have naively thought in the case the dust is immediately 
released into the ISM a short time after supernova explosions of 
massive, short-lived stars which dominate the SFR. These two 
peaks are separated by ~ 2.5 Gyr, which is a long period if 
one considers that the dust production rate peak is below 1 Gyr 
for SNII, bu t it is at 3-4 Gyr for intermediate-mass stars (see 
iDwekl 19981) . Nevertheless the dust reaches a sort of plateau from 
z ~ 1 .5 (with a m aximum at z ~ 1) , i.e. 1 Gyr after the peak of the 
SFRD. Recently, iFukugital (1201 ll) suggested that the dust must 
survive on much longer time scales than what has been previ- 
ously thought and that half of the dust could be produced by SNII 
and the other half by intermediate-mass (1-8 solar masses), long- 
lived stars. If we assume that the SNII dust production peaks 
very shortly after the SFRD peak, then the dust peak that we 
observe at z ~ 1 is likely due to intermediate-mass, long-lived 
stars producing their peak of dust on a delayed time. In partic- 
ular, the AGB stars release dust through intense mas s loss, and 
most efficiently at the very end stages of evolution (Gall et al. 
1201 ll) . The peak of attenuation is higher than SNII simply be- 
cause these stars are much more numerous than very massive, 
short-lived OB stars, assuming a universal IMF. It would explain 
also that at z > 2, the low level of dust attenuation is mainly 
due to dust produced by SNII, while at z < 2, it is resulting 
from the combination of dust from SNII and intermediate-mass 
stars. Our findings, combined with the above-mentioned times 
scales, may imply that dust is not only quickly formed (thanks to 
SNII ejecta and remnants), but also quickly destroyed, because 
the peak of attenuation is only ~ 2.5Gyr after the SFRD peak, 
a shorter time-scale than the typical one for dust production by 
intermediate-mass stars (3-4 Gyr). Surely, like in the case of the 
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SFRD evolution, our findings about the general evolution of the 
FUV-band dust attenuation in such a broad cosmic epoch (~ 12 
Gyr) will constitute an important reference for future models. 
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Appendix A: The weighting schemes 

We refer the reader to iLe Fevre et al 1 (I2005bl) and Le Fevre et 
al. (in prep.) for details on the quality flags assigned to redshift 
measurements in the VVDS Deep and Ultra-Deep surveys, re- 
spectively. We detail below the weighting schemes to be applied 
to our statistical analyses. They account for the selection func- 
tion of the photometric sources targeted to acquire their spectrum 
(the Target Sampling Rate, TSR), and for the success to measure 
a reliable redshift from the spectrum (the Spectroscopic Success 
Rate, SSR). They are derived for the Deep survey, the Ultra- 
Deep survey, and for the merged Deep+Ultra-Deep sample. 

To derive the weighting schemes, we made use of the pho- 
tometr ic redshifts, zphot, computed as described in lllbert et ail 
(120061) . but using the more recent T0005 release of CFHTLS 
photometric data (u*,g',r',i',z r ) and the latest near-infrared 
photometric data available from WIRDS (J, H and K s , 
iBielbv et al.ll201ll> . 

We emphasise that we correct the Deep and Ultra-Deep dis- 
tributions of photometric redshifts, n{zphot), for their failure rate 
computing the ratio between the spectroscopic redshifts mea- 
sured with a > 97% confidence level (flags 3, 4), «(z/=3,4) an d 
their corresponding photometric redshifts, n(zphotf=x4)- 

A.1. Weights for the Deep survey 

The Deep survey TSR^ is defined as Nf arget /N^ hot , where N d phot 
is the number of photometric sources in the Deep survey with 
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17.5 < Iab - 24.0, over ~2200 arcmin 2 of sky area, and Nf a t 
is the number of photometric sources targeted for spectroscopic 
observations. The TSR^ depends on the projection of the angular 
size of each objec t on the x-axis of the image ('x-radius', see also 
lllbert et alj|2005t) . The weight associated to TSR,/ is wt„„(r) — 
1 /TSRd(r), where r is the x-radius defined above. 

The Deep survey SSRj is defined as N d pec INf arget , where 
N d c is the number of targets with a reliable spectroscopic red- 
shift measurement in the Deep survey. The SSRj is a functio n 
of both selection magnitude and redshift (see lllbert et al.ll2005l) . 
As described in Le Fevre et al. (in prep.), a fraction of low con- 
fidence level spectroscopic redshifts (flags 0, 1, 2) at z > 1-4 
from the Deep survey have been observed again, leading to a 
100% confidence level in the redshift measurement. Comparing 
the old redshift distribution n{z'J°y2 >ved ) to the new one for this 
subsample, we have remodulated the full n(z/ =1,2) partially us- 
ing the n(zphot) for z < 1.4 sources. We have also remodulated 
the single emission-line redshifts (flags 9) n(zf=9) according to 
its n(zphotf = g). This gives: 



0.4 F 
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Fig.A.l. Target Sampling Rate (TSR) for the merged 
Deep+Ultra-Deep catalogue as a function of the /-band appar- 
ent magnitude. Given the different flux limit of the two surveys 
(see text), this TSR is obtained using Deep data at Iab < 23, both 
Deep and Ultra-Deep data at 23 < Iab ^ 24, and Ultra-Deep data 
at Iab > 24. The last point is low because the Ultra-Deep survey 
TSR is much lower than the Deep survey TSR. 



N d spec (m,z) = N f=ia (m,z)xn(m,zjf l T ved )/n(m,z f=1 ,2) + 
+ Nf = g(m,z) X n(m,zphotf = g)/n(m,Zf=9) + 
+ Nf = 3,4(m,z). 

The weight associated to SSR</ is vf' SSR (m, z) = 1/S S Rd(m,z). 

In summary, for the Deep survey, we have Nf pec c Nf arget c 
N d phot and we apply the weight Wd(r,m,z) = w^ssM w d SSR (m,z) 
to each galaxy in Nf pec . 

A.2. Weights for the Ultra-Deep survey 

The Ultra-Deep survey covers a 576 arcmin 2 sky area embed- 
ded in the Deep survey area. It is purely flux limited at 23.00 < 
i' AB < 24.75 with N p d wt photometric sources. Several spectro- 
scopic sources with 23 < Iab ^ 24 had already been observed 
in the VVDS Deep survey, when we started the Ultra-Deep 
observations. By excluding them, we obtained a reduced cata- 
logue with N p ^'™' sources available for the Ultra-Deep spectro- 
scopic target selection. The Ultra-Deep Photometric Sampling 
Rate (PSR„ rf ) is defined as N p ph r o e "' IN" d wt . The PSR„ d depends on 
the /-band apparent magnitude. The weight associated to PSR„</ 
is w ud SR {m) = l/PSRudim). 

The Ultra-Deep survey TSR lld is defined as N'/^JN^"' , 

where N" d rget is the number of photometric sources targeted for 
Ultra-Deep spectroscopic observations. The TSR„</ does not de- 
pend on any parameter, it is a constant value (6.5%). The weight 
associated to TSR U£ / is w l j SR = l/TSR u d- 

The Ultra-Deep SSR ud is defined as Nf p JN l ^ rge „ where 
N" pec is the number of targets with a reliable spectroscopic 
redshift measurement in the Ultra-Deep survey. Here N spec — 
N/= 1.5,2,3,4,9, where flags 1.5 is a flag 1 with a photometric red- 
shift in agreement with the spectroscopic redshift (Le Fevre et 
al., in prep.). The SSR„ is a function of both selection mag- 
nitude and redshift. The associated weight is w ud SR (m,z) = 



l/SSRud(m,z). 



1 target 



In summary, for the Ultra-Deep survey, we have N" pec c 



c N pZT c N phot and we aPPty the wei 8 ht w ud (m, z) 



w% R (m) w% R w SSR 



ud „,ud t m ^ z ) tQ g^jj g a J aX y m N" pec . 



A.3. Weights for the combined Deep and Ultra-Deep surveys 

We combine the two surveys and we derive an adapted weighting 
scheme to take advantage of both (i) the large magnitude range 
covered by the Deep survey (17.5 < Iab ^ 24.0) over a large sky 
area, and (ii) the depth reached by the Ultra-Deep survey (0.75 
mag deeper than the Deep) over a smaller embedded sky area. 

As a first step, we need to merge the two catalogues of galax- 
ies in accounting for their different depth, especially for their dif- 
ferent LF bias limit (see Sect. 13.11 ). Since the latter depends on 
the studied redshift bin, we merge the catalogues in each redshift 
bin in which we explore the LF, as follows. For the Ultra-Deep 
survey, we compute the LF bias limits with the method described 
in Sect. 13.11 For the Deep survey, we adopted more conservative 
limits with respect to those obtained with the above mentioned 
method, because of the magnitude range 23 < Iab ^ 24 com- 
mon with the Ultra-Deep survey. In the Deep and Ultra-Deep 
samples we keep only galaxies brighter than the respective bias 
limit, and we merge the two sub-samples into a single catalogue. 
From now on, the merged catalogue is considered with flux lim- 
its within 17.50 < Iab ^ 24.75, over the Deep survey effective 
sky area. 

As a second step, we need to slightly modify the TSR 
weighting scheme for the merged catalogue since the Ultra- 
Deep area is embedded in the Deep area. The Deep+Ultra-Deep 

TSR, +1 „ is (Ni rget + N» d rge ,)IN f p f ot , where iv£ is the full photo- 

metric catalogue at 17.50 < Iab ^ 24.75 over ~ 2200 arcmin 2 . 
The TSR u+u d depends on the /-band apparent magnitude (see 
Fig. IA.lt . At 17.5 < I A b < 23.0, it corresponds to TSR rf when 
we compute it as a function of magnitude instead of angular size, 
while at 24.00 < Iab ^ 24.75, it corresponds to the constant 
TSR„ f /. At 23 < Iab ^ 24, it corresponds to the combination of 
both TSRs, i.e. at the numerator we have all targets of both the 
Deep and Ultra-Deep surveys. The weight associated to TSR f /+,„/ 
is Wj + s ' R (m) = l/TSRd+ud(tn). We applied this weight according 
to the /-band apparent magnitude to all Deep survey galaxies 
and to Ultra-Deep galaxies which are brighter than the LF bias 
limits of the Deep survey. To the Ultra-Deep galaxies which are 
fainter than the Deep LF bias limits we applied the w'j + s " R shown 
for/ Afl >24. 
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Fig.A.2. Rest-frame FUV l/Vmax LF data points for the Deep 
survey (cyan empty squares) and the Ultra-Deep survey (blue 
filled triangles), from z = 0.05 to z — 4.50. There is no Ultra- 
Deep data in the lowest redshift bin. In each panel, the FUV 
absolute magnitude LF bias limit (see Sect. 13. Il l is shown as a 
vertical cyan long-dashed line for the Deep survey and as a blue 
short-dashed line for the Ultra-Deep survey. The number quoted 
in parenthesis is the total number of used galaxies. We overplot 
the STY determinations of the merged Deep+Ultra-Deep sam- 
ple, as illustrated in Fig. Q] and tabulated in Table IBTTT i: (i) with 
free Schechter parameters (solid lines) and its associated error 
(shaded area), (ii) in setting the faint-end slope at z > 1.7 (dashed 
lines), and (iii) in setting M* FUV at z < 0.2 (dashed line) and re- 
ported as a reference in the z > 0.2 panels (dotted lines). 



In summary, for the merged Deep+Ultra-Deep catalogue, we 
apply the weight w™ erged (m,z) = w d T + s ud (m) w d SR (m,z) to N d spec 
and the weight w™ rged (m,z) = wp s ^(m) w ltd SR (m,z) to N$ ec . 

A.4. The FUV-band LFs using the different weighting 
schemes 

To verify the robustness of our weights, we separately derive 
the rest-frame FUV LF of the Deep and the Ultra-Deep sur- 
veys, using their independent weighting schemes, as described 
in Appendixes I A. 1 1 and IA.2I Figure I7v2l shows the l/V maA data 
points up to the respective LF bias limit (see Sect. 13. Il l for the 
two surveys, together with the STY LFs of the combined Deep 
and Ultra-Deep surveys (see Sect. 13.21 and Fig. [TJ, using the 
merged weighting scheme, as described in Appendix IA.3I On 
one side, the LF data points of the two surveys fully agree from 
z = 0.2 to z — 3.5 in their common range of luminosities. This is 
an a-posteriori check of the reliability of our individual weight- 
ing schemes. On the other side, the STY LFs estimated with 
the merged sample perfectly overlap the 1/V max data points es- 
timated with the two individual surveys. This confirms the relia- 
bility of the weighting scheme applied to our merged sample of 
individual surveys with different depth. 



Fig. B.l. Rest-frame FUV l/V max LF estimates from z = 0.05 to 
z = 4.5 derived with the VVDS zphot catalogue up to Iab = 26 
(red open diamonds). They are plotted over a grey-scale ver- 
sion of Fig. [H that is the rest-frame FUV LF of the VVDS 
Deep+Ultra-Deep spectroscopic sample up to Iab = 24.75. 



Appendix B: Verifying the robustness of our 
rest-frame FUV-band LF 

In Sect. [3] we have presented our VVDS rest-frame FUV-band 
LF estimates from z = 0.05 up to z = 4.5 (see Fig. [TJ, and we 
have tabulated the corresponding best Schechter parameters in 
TableQ] Figure|2]summarises in one panel our best LF fits. In this 
Appendix, we give further details and tests we have performed 
to ensure the robustness of our best rest-frame FUV LF, in par- 
ticular the flatness of the LF faint-end slope at z S> 1 and the low 
normalisation of the LF at z ^ 3. Table lETTI gives our extensive 
computations of the Schechter parameters with the STY method, 
leaving the three Schechter parameters totally free. Obviously, 
we cannot constrain M* at z < 0.2 nor a at z > 1.7. We fixed 
these two parameters as described in Sect. 13.21 and we also re- 
port the results in Table IB. II The entries labelled with a black 
dot are our final choices taken for our studies and also reported 
in Table [lj 

As shown in Fig. [2] our FUV faint-end slope is quite flat 
(or ^ -1) at z < 1.7, that is up to the highest redshift where we 
can constrain it. The compilation of the a values found in the 
literature (see Table IB.2"l i shows very scattered values at < z < 
6. Nevertheless, it is a notable point that a has been generally 
found steeper (when estimated and not fixed) than our value. 
To strengthen the validity of our results against possible biases, 
we report below our LFs computed using photometric redshifts, 
using only starburst galaxies, and computed in the rest-frame 
NUV-band. 

B.1. The rest-frame FUV LF with zphot 

We derive the VVDS rest-frame FUV l/V„ wx estimates using 
very deep photometric redshifts (zphot, see Appendix [A] for de- 
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Table B.l. Schechter parameters (M*, a, (f>*) and total Luminosity Density LD for the rest-frame FUV- and NUV-band galaxy 
luminosity function for the VVDS Deep+Ultra-Deep sample with the cosmology (Q m , Q.^, h) = (0.3, 0.7, 0.7). 











FUV 








NUV 






Redshift range " 


M* 


a 


<t>* 


lg(LD Ilt .) * 


M* 




0* 


lg(LD„ c ) " 






(AB mag) 




(IQ- 3 Mpc~ 3 ) 


W/Hz/Mpc 3 


(AB mag) 




(10~ 3 Mpc" 3 ) 


W/Hz/Mpc 3 




Free Schechter parameters 




0.05 < z < 0.2 
0.2 < z < 0.4 
0.4 < z < 0.6 


-18.6itf 
-18.4* 


_! j +o.o6 
i 17+8:8?? 

-107*^ 


5.36+M6 
6.91+1:02 

6 60*^ 

"•""-0.86 


18.86-- 

18.87+; 

18.85+; 


-18.5+; - 
-18.3* - 
-18.7* - 


1 ns+o- 06 

LUS -0.0fi 

1 02+"- !i + 

L08 -i 


6 35+' M 
10.16!j 
6.S2+ 4 
9 27+' " 1 
9 48+ o: ^ 7 

7 -^°-0.88 

6 57+0- 86 

"•-"-0.89 

3 49+0.58 
q 7^3+1 :^9 

J - /J -1.67 
1 09+0.37 
1 -- ,z '-0.47 


18.90!-' 
18 97+ 8 - 03 

18 99+° .02 

19 12+"- !i2 

-Q-P2 




0.6 < z < 0.8 


-18.3* 


-0.90*|| 


Q so +0.99 
-Q.SS 


18.93* 


-18.8* - 


0-95* 




0.8 <z< 1.0 
1.0 < z < 1.2 
1.2 <z< 1.7 
1.7<z<2.5 
2.5 <z<3.5 


-18.7*! 
-19 0+ 02 

17 '"-0.2 

-19 6 +02 

ly -"-0.2 

-20.7+ , 4 
-20.8*3 


-0-85*1 
-0 91+ 016 

"• 71 -0.16 

-1 09+° 23 
-1 65+ 055 
-0 63+ 057 

"•"-'-0.53 


9 01 +!i?4 

7 -" x -0.96 
7 43+1-08 
'•^-'-1.15 

4 10+ 077 

9 9Q + 1.85 

• J -l£6 
i 79+O.19 
1,/z, -0.36 


19 04+0 02 

17 ' u -0.01 
1 Q 1 9+0.04 
-0.03 

19 13+ 012 

17,1 -0.07 

19 70+ 3 25 

"■'"-0.36 

19 17+ 025 

iy l '-0.10 


-19.2*! - 

_19 6+oi 
17 -"-o.i 

-20.2+01 - 

-20.7+o :3 - 

-21 4+ 03 
i1- -0.3 


0.81*| 
0.88!» : | 3 
1 05 + 017 

I -" J -0.16 

1 16+ - 47 

I - I "-0.45 

1 15+ 042 
1 - 1J -0.40 


1 Q 97+0.01 
17 -^'_0.01 

19 30+ 02 
19 28+ 06 
19 56+° 52 
19 37+ - 38 

i7 -'"-0.13 




3.5 <z<4.5 




















M* constrained , a and <p* free 


• 


0.05 < z < 0.2 


-18.12 


-1 05+ U04 

1 -" J -0 04 


7 00+"" 44 
'•""-0 44 


18.76+0- 


















a constrained 


M* and 0* free 








• 
• 
• 


1.7<z<2.5 
2.5<z<3.5 
3.5 < z < 4.5 
1.7<z<2.5 
2.5 < z < 3.5 
3.5 < z < 4.5 
1.7<z<2.5 
2.5 < z < 3.5 


-20.4—j 
-21 4 +o:i 

97 9+O.2 

-8? 

-20 3 +iiT 
-Q.J 

-2i.ii 

-21.8*! 
-20.8+ ,; 
-21 6 +!i:T 


-1.30 
-1.50 
-1.73 
-1.10 
-1.10 
-1.10 
-1.73 
-1.73 


Q ^7+0.24 

0.86!;| 

0.1 1+ 00 ' 
3 94+8:24 

i 97+8:^? 

• '-Q.Q7 

22+ - 02 

"•^-0.02 

1 95+014 
60*^ 

"•""-0 03 


19 46+o oJ 
i^.tu _ 003 

19 40+o .02 
19 10+o .06 
19 41+0 05 
-QQ2 
19.245 

18.73+™ 3 
19 79+0.05 

19 62+^ : ^ 











Notes. (fl) The entries labelled with a filled dot are the chosen Schechter parameters taken for the rest-frame FUV LF used throughout the paper, 
and reported in Table[TJ <fc) The quoted errors are those induced by the STY LF fit; see T ableQJfor errors including also other sources of uncertainty. 
W y[* j s se [ at z < 0.2 assuming the local rest-frame FUV value determined bv lWvder et al.ll2005i) . <rf) a is set with different values at z > 1.7. First, 
following an evolution with redshift as described in Sect. l3.2l (g = -1.3, -1.5, -1.73 at 1.7 < z < 2.5, 2.5 < z < 3.5, 3.5 < z < 4.5, respectively); 
secondly, set as two likely extreme non-evolving values (a = -1.1 and a = -1.73). We remind that a = -1 .1 corresponds t o wha t we find at 
1.2 < z < 1.7 (where we can still determine a with a reasonable error), and a = -1.73 is the value found bv lBouwens et alj d2007f) at z = 3.8. 
Contrarily to the rest-frame FUV LF, we were not interested in setting any Schechter parameters for the rest-frame NUV LF. 



tails). The zphot catalogue is complete up to Iab = 26 and it 
consists in ~ 43000 photometric sources over the VVDS Deep 
area. In this case, the weighting scheme is obviously not used. 
The bias limit has been computed in each z bin as explained in 
Sec. 13. II and at any z it is clearly fainter than the bias limit in the 
VVDS Deep+Ultra-Deep spectroscopic sample. 



Results are shown in Fig. IB. II and globally overlap our 
best rest-frame FUV-band LF determinations using the VVDS 
Deep+Ultra-Deep spectroscopic sample up to Iab = 24.75. 
Thus, our spectroscopic sample (brighter 1.25 mag in Iab than 
the zphot sample) does not significantly miss any faint galax- 
ies, which strengthens the reliability of our flat faint-end slope. 
One small discrepancy occurs at z < 0.2, where the photomet- 
ric 1/Vma.i points have an higher normalisation with a slightly 
steeper slope. It is due to the well known degeneracy in the 
computation of photometric redshifts, i.e., some objects end up 
wit h a wrong too low (zphot < 0.3) photometric redshift (see 
e.g. Illbert et al.ll2009h . Furthermore, it confirms once more that 
our weighting scheme is adequate. This strengthens also the 
reliability of our low (p* value at 3.5 < z < 4.5, which is 
neither caused by the cut at Iab = 24.75 nor by an incom- 
pl ete weighting scheme. W e refer the reader to Sec. 13.3.21 and 
to iMcCracken et al.l (120031) for a discussion about photometric 
completeness with respect to /^B-band magnitude and surface 
brightness. 



B.2. The rest-frame FUV LF of starburst galaxies 



We are interested to check whether the flatness of our slope at 
z < 1 .7 could be due to a very low density of intrinsically reddish 
galaxies, which could have hidden, in the total sample, very high 
densities of faint intrinsically blue galaxies. Indeed an /-band se- 
lection does not select only starbursting g alaxies, in contrast w ith 
f/V-band selected samples like GALEX (lArnouts et al.ll2005l) . 

We have classified our galaxies of the VVDS Deep+Ultra- 
Deep spectroscopic redshift sample according the VVDS 
schem e of four photometric types, as described in IZucca etall 
Here, we consider the 'Type 4' galaxies, correspond- 
ing to galaxies for which the best template fit over the 
u*g'r'i'z'JHK s broad bands is among the bluest galaxy tem- 
plates, i.e., a starburst or irregular galaxy template. The Type 4 
rest-frame FUV-band 1 / V nm x LF estimates are shown in Fig. IB. 21 
At z < 1 -7, their normalisation is obviously lower with respect to 
the total sample, because the fraction of Type 4 galaxies in our 
sample goes from ~ 40% at z ~ 0.1 up to ~ 85% at z ~ 4 (when 
considering galaxies brighter than the bias limit). The faint-end 
LF slope for Type 4 galaxies is very similar to the one found 
for the total sample. Our flat slope is obviously not due to very 
low densities of faint intrinsically red galaxies, which could have 
overcompensate very high densities of faint intrinsically blue 
galaxies. Thus, it reinforces the flatness of our slope. 



19 



O. Cucciati, L. Tresse, O. Ilbert et al.: The CSFRD since 12 Gyr traced by the VVDS 




-2 
-4 

-6 
-8 

-2 
-4 

_6 

Q. 

? -2 
W -4 
1 -6 
s. -8 

o _ 2 
-4 
-6 
-8 

-2 
-4 
-6 



0.05-0.20 — T27BT 




is e O o o Q 



0.40-0.60 — [TTOTJ 




0.80-1.00 — (1611) 




1 .20- 1 .70 i ( 976) 




2.5d-3!50 (316) 




0.20-0.40 — (1219} 



0.60-0:S0 — 




1.00-1:20 — T9407 




1.70-2:50 — ( 257) 



3.50-4:60 — ( 130) 



-25 



-20 



-15 



-25 



-20 



-15 



Fig.B.2. Rest-frame FUV l/V,„ ax LF estimates from z = 0.05 
to z — 4.5 for our intrinsically bluest galaxies (red open cir- 
cles, VVDS Type 4, see text) in the VVDS Deep+Ultra-Deep 
spectroscopic sample. They are plotted over a grey-scale ver- 
sion of Fig. [H that is the rest-frame FUV LF of the total VVDS 
Deep+Ultra-Deep spectroscopic sample. 



B.3. The rest-frame NUV LF 

The last point to test is the fact that our rest-frame FUV mag- 
nitudes are based on template extrapolations at z < 1, which 
could have an impact on the faint-end LF slope. Our bluest pho- 
tometric information is in the m*-381 CFHTLS filter, and thus 
the rest-frame FUV emission is directly observed from z > 1 . 
Nevertheless, the rest-frame NUV emission is directly observed 
at z > 0.4. Thus we report our results from the rest-frame 
NUV-band LF of the Deep+Ultra-Deep spectroscopic sample. 
We list the Schechter parameters derived with the STY method, 
in Table IB. fl and plot the LFs in Fig. IB. 31 We note that we were 
not interested in trying to fix M* at low z, nor to fix a at very 
high z. These estimates enable us to measure the faint-end slope 
at 0.4 < z < 1, that is where the FUV intrinsic luminosity is 
extrapolated but not the NUV intrinsic luminosity. Again, we do 
not find a steeper faint-end slope in the rest-frame NUV-band 
LFs. As our aim is to derive the SFRD from the rest-frame ul- 
traviolet c ontinuum spannin g from FUV to the NUV using the 
relation of lKennicutj (1 1 998b . we can affirm that the extrapolation 
of the FUV-band luminosities does not change our SFRD results 
since the faint-end slope stays flat. 

Appendix C: FUV number counts and luminosity 
functions at z ~ 3 

At 2.7 < z < 3.4. lLe Fevre et al.l (l2005ah and Le Fevre et al. (in 
prep.) show that the VVDS presents number counts per 0.5 mag 
interval of Iab apparent magnitude per unit surface area (i.e., 
dN/0 . 5 mag/ arcmin 2 ) at least 2 times larger than those quoted by 
ISteidel et al.l ( 1 19991) within the same redshift range. This is par- 
ticularly evident for bright galaxies (Lib < 23.5). In contrast, the 



Fig.B.3. Rest-frame NUV-band LFs from z = 0.05 to z = 4.5 
of the VVDS Deep+Ultra-Deep sample, in 10 redshift bins as 
indicated in each panel. Red circles represent the 1/V max data 
points up to the LF bias limit and the corresponding number of 
galaxies are given in parenthesis. The vertical cyan long-dashed 
and blue short-dashed lines correspond to the LF bias limits for 
the Deep and Ultra-Deep surveys, respectively. The black solid 
curve and the associated shaded area is the STY LF estimate 
assuming free Schechter parameters and its associated error. The 
dot-dashed curve is the same in each panel, and it corresponds to 
the solid curve in the 0.2 < z < 0.4 redshift bin, as reference. All 
Schechter parameters are given in Table IB. ll For the NUV-band, 
we were not interested in fixing either M* nor a. We note that in 
the last redshift bin the STY fit is unconstrained, so there is no 
solid curve overplotted. 



VVDS re st-frame FUV LF is only ~ 50% higher than the one 
found by Steidel et al. ( 1999), as illustrated in the left panel of 
our Fig. IC.ll and this holds only for Mfuv < -22.5. Here, we 
demonstrate that these two results are not in conflict. We remark 
that throughout the paper we compute our VVDS FUV absolute 
magnitudes using a filter centred at 1500 A, while ISteidel et al.l 
(fl999l) use a FUV-band filter centred at 1700 A. To compare our 
results with their work, in this section we comput e our VVDS 
FUV a bsolute magnitudes using the same filter as ISteidel et al.l 
( fl999l) . even though our VVDS FUV-150 and VVDS-1700 LFs 

are ver y similar. 

As ISteidel et alJ dl999h begin their analysis from galaxy 
number counts as a function of % apparent magnitude, we over- 
plot the same for the VVDS Deep and Ultra-Deep surveys in 
the right panel of Fig. IC.ll We remark that in this Figure we 
plot the number of galaxies per unit magnitude and per unit vol- 
ume (i.e., dN/mag/Mpc 3 ), while lSteidel et all d 19991) cite in their 
Table 3 the surface density number counts for 0.5 mag intervals 
(dN/0.5mag/arcmin 2 ). To convert their values in dN/mag/Mpc 3 
units, we multiplied them x2 to obtain the correct magnitude in- 
terval, and we divided for the effective volume of each arcmin 2 , 
quoted in their Table 3 (with the same cosmology that we adopt, 
Cl m = 0.3, Qa = 0.7). Our number counts are weighted as de- 
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Table B.2. Non exhaustive list of the faint-end slope values, a, of the rest-frame FUV-band luminosity functions at < z < 6, as 
found in the literature. Values are sorted according to the redshift. The effective wavelength of the rest-frame band at which the LF 
was determined is also given. 
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Notes. (n> Mean or median redshift as quoted in each work, or, if not specified, centre of the studied redshift interval. {b) Faint-end slope of the UV 
LF, with error bars when it was estimated, with n o error bars whe n it was fixed. ic) Faint-end slope of the nearby FUV LF estimated after setting 
M* FUV . (rf) Value retrieved from the compilation of Hopkins (2 0041) . 



scribed in Appendix l A. 1 l and lA. 21 and they are realized as a func- 
tion of msT = (r + 0/2 apparent magnitude, i.e., the mean value 
of the r- and /-CFHTLS br oadband filters, beca use it mimics the 
•R-broadband filter used bv lSteidel et al.l d!999l) . We observe that 
the VV DS counts are clearly higher than those from lSteidel et al.l 

(fl999h . at least at the brightest magn itudes. 

Next, from their number counts, ISteidel et al.l (1 19991) derive 
their rest-frame FUV LF (shown in the left panel of Fig. IC.lt . 
translating apparent magnitudes into absolute magnitudes with 
the following rigid shift: 

M FUV = m K - DM(z med ) + 2.5 log(z med + 1), (C. 1) 

where DM(z me d) is the distance modulus defined as 5 x 
log(D£(z me rf)/10pc), with Di{z me d) the luminosity distance at 
Zmed, and Zmed the median redshift of the studied sample. For their 
sample at 2.7 < z < 3.4 {z me d = 3.04), it gives: 

M FUV = m K -45.55. (C.2) 

For the VVDS Deep+Ultra-Deep sample, at 2.7 < z < 3.4 
we have z me d = 2.92. With the same method, starting from our 



observed number counts (see right panel of Fig. IC.lt . we obtain 
the rest-frame FUV LF data points, as illustrated in the left panel 
of Fig. IC.ll and we fitted them with a Schechter function. The 
result is that our c/JV-derived LF is higher at the bright end than 
our direct 1 / V max LF data points, as derived in Sect. 3.2, and here 
fitted also with a Schechter function. We conclude that our c/JV- 
derived LF presents a different shape than our 1 / V mflx -derived 
LF. 

Below we verify that, as expected, it is not correct to use a 
simple, rigid shift to transform our apparent magnitudes to ab- 
solute ones, as done with Eq. IC.ll to obtain correct c/JV-derived 
LFs. The main reasons are the following. 

1) The n(z) shape. In the range 2.7 < z < 3.4, the redshift distri- 
bution n(z) of an /-band selected survey as observed by the 
VVDS decreases as a function of z as shown in Le Fevre 
et al. (in prep.), in particular in a faster way for brighter 
galaxies. In contrast, due to the colour selec tion function 
used t o find Lyman-break galaxies, the n(z) in ISteidel et al.l 
( 1 19991) peaks around its median value within the same red- 
shift range. The shapes of the different n(z) are shown in the 
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Fig. C.l. Left panel. Comoving luminosity densities as a function of FUV (1700 A) absolute luminosities at 2.7 < z < 3.4: 1 /V max 
estimates of our VVDS Deep+Ultra-Deep sample fitted with a Schechter function (red circles and dashed curve, respectively); our 
nisT = (r + 0/2 number counts per volume density for the Deep and Ultra-Deep surveys converted into FUV absolute magnitudes 
with a rigid shift (see text), and fitted with a Schechter function ( blue diamonds and triangles, and solid line, respectively); the 
nisT = H number counts per volume density of ISteidel et al.l (1 19991) converted into FUV absolute magnitudes with their rigid shift 
(see text), and fitted with a Schechter function by them (orange squares and solid line, respectively). Right panel. Number counts 
per unit volume of galaxies as a function of msr apparent magnitudes at 2.7 < z < 3.4: our number counts from the VVDS Deep 
and Ultra-Deep surveys (blue diamonds and triangles, respectively); our VVDS number c ounts derived from o ur VVDS 1 /V max LF 
(see right panel) shifted rigidly in magnitude (red circles, see text); number counts from lSteidel et al.l (1 19991) (orange squares; see 
text for conversion details). The three curves correspond to the sa me curves in t he left panel, shifted rigidly (see text). Note that the 
CFHTLS (r + i)/2 photometric system mimics the fi one used bv lSteidel et al.ldl999l) . 



inset of Fig. lC2l If the n(z) is skewed, it is crucial to use the 
correct redshift distribution in Eq. IC.ll and not simply the 
median redshift value. 
2) The K- and colour corrections. ISteidel et alJ (1 19991) restrict 
the (/f+colour) term to 2.5log(l + z), as in Eq. IC.ll because 
they work with the observed 'R-band, corresponding to the 
FUV light emitted at z ~ 3. Nevertheless, this is correct for 
z ~ 3.0 galaxies only. 

We made a simple simulation to illustrate these two aspects. 
We started with a set of 3 x 10 6 galaxies distributed in FUV 
absolute magnitude like the Schechter fit of our l/V, !M V -derived 
FUV LF. We want to know whether, with a selection function 
as the VVDS, we recover our observed galaxy number counts as 
a function of nisr apparent magnitude at 2.7 < z < 3.4. Using 
Eq. IC.ll we operate a rigid shift to transform our FUV abso- 
lute magnitudes into ntsr apparent magnitudes, and to obtain 
a LF-derived number counts as illustrated in the right panel of 
Fig. IC.ll We call this LF-derived distribution n(msT, LF). We 
see that these LF-derived number counts are not in agreement 
with our observed number counts which are at least twice higher 
at nisr < 23. We illustrate in Fig. lC.2l the ratio of our observed 
number counts to our n(msr, LF) as a function of m^j. 

We will now test the n(z) and the (/f+colour) corrections, 
respectively. We fit separately the n(z) of VVDS galaxies with 
nisr < 23 and with nisr > 23, to account for the steeper n(z) 
decrease for bright magnitudes (see inset in Fig. IC.21 >. Given the 
monotonic (even if scattered) relation between observed and ab- 
solute magnitudes, ntsr ~ 23 corresponds to Mfuy ~ -21.1. 
Thus, for each previously simulated galaxy, we assign a redshift 
in respecting the different n(z) distributions for galaxies brighter 
or fainter than Mfuv = -21.1, that is we do not anymore as- 



sume a single median redshift. We now compute the observed 
msr magnitudes of our simulated sample using Eq. IC.ll where 
we substitute z me d with the particular redshift that has been as- 
signed to each galaxy. We obtain a msr distribution that we 
call n(msT,z). The ratio of our observed number counts to the 
n(msT,z) is shown in Fig. IC.2I The n(msT,z) distribution is now 
closer to our observed number counts than the n(msr,LF) dis- 
tribution. 

As a further step, we account for the (K+colom) terms for 
each previously simulated galaxy, that is: 

M FUV = m ST - DM(z) - K ST (z) - (m ST - m FUV ) (SED ' z=0 \ (C.3) 

where Kst(z) is the ^-correction at z and (nisr - mpuv) is 
the colour term at z — 0, both derived from the best SED fitting 
template. We have modelled the best fitting template distribution 
in the VVDS sample, which gives us the VVDS distributions of 
Kst(z) and (nisr - m.fuv) ( - SED ' z={y> terms. Also in this case we 
assumed the two n(z) distributions (for faint and bright galaxies) 
described above. Using the specific distribution of all the terms 
in Eq. IC.3I we computed a new set of msr starting from our 
sample of absolute magnitudes, and obtained their distribution 
{n(msT, z, K, col)). We show in Fig. lC.2l the ratio of our observed 
number counts to the n(msT,z, K,col) distribution. We observe 
that the modelisation of the A'+colour term further improves the 
match with the observed number counts. Still, the ratio is not 
equal to unity, which demonstrates that the reality is more com- 
plex than a simple simulated recipe. 

To mimic the VVDS observed data, but still with a very sim- 
ple test, we have fitted with a Gaussian function the distribu- 
tion of AM = Mfuv - msr in the VVDS sample. Subtracting 
from each simulated Mfuv a value of AM extracted randomly 
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Fig. C.2. As a function of msr, ratio of the directly observed 
VVDS FUV number counts at 2.7 < z < 3.4 to the number 
counts derived in various ways from the FUV luminosity func- 
tion in the same z range. The black dotted line (ratio=l) is the 
ideal case where both counts are identical. The black solid line 
represents the ratio of the observed VVDS to the VVDS LF de- 
rived counts using Eq. IC.ll nfmsr, LF) (i.e. blue solid line to the 
red dashed line of Fig. IC.ll right panel). The other lines repre- 
sent the ratio of the observed VVDS number counts (the above 
mentioned blue line) to the following various simulated distribu- 
tions: n(msT,z) (blue long-dash curve), n{nisT,z,K,col) (green 
short-dash curve) and n{msT, AM) (red dot-dash curve). See text 
for details about these simulated distributions. In the inset, we 
plot th e global shape of the redshift distributions in lSteidel et al.l 
d 19991) (dashed line) and in our work, separately for faint and 
bright galaxies (thin and thick solid lines, respectively). The nor- 
malisation on the y-axis is consistent between the two solid lines, 
while the normalisation of the dashed line has been arbitrarily 
chosen to better illustrate the different shapes. 



to|Le Fevre et aT] (12005 al) and Le Fevre et al. (in prep.) for more 
details on number counts and sample selections. 
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from the distribution, we obtained a new nisr distribution, that 
we called n(msr, AM). It is clear that this last distribution (see 
Fig. IC.2t mimics better than the others the real observed number 
counts distribution at bright magnitudes, but then the match for 
fainter galaxies is slightly worse. 

In summary, this simple exercise confirms that it is very dan- 
gerous to transform number counts within a given redshift range 
to a luminosity function in the case of a skewed n(z) within 
the redshift interval considered, and/or in the case of dissimi- 
lar n{z) for different galaxy populations. The fact that the n{z) of 
colour-selected samples is very different from the n{z) of mag- 
nitude limited samples, implies that number counts when us- 
ing colour-selected samples are not representative of the number 
counts for the complete galaxy population. We refer the reader 
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